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ABSTRACT 

We describe the physical and orbital properties of C/2011 W3. After surviving the perihelion pas- 
sage, the comet was observed to undergo major physical changes. The permanent loss of the nuclear 
condensation and the formation of a narrow spine tail was observed first at Malargue, Argentina, on 
December 20 and then systematically at Siding Spring, Australia. The process of disintegration culmi- 
nated with an outburst (terminal fragmentation event) on December 17.6 UT. The postperihelion tail, 
observed for ^3 months, was the product of activity over <2 days. Because of the delayed response to 
the hostile environment in the immediate proximity of the Sun, the nucleus' breakup and crumbling 
was probably caused by thermal stress due to the penetration of intense heat pulse deep into the nu- 
cleus' interior after perihelion. The same mechanism may be responsible for cascading fragmentation 
of sungrazers at large heliocentric distances. The observed behavior is at odds with the rubble-pile 
model, since the residual mass of the nucleus after perihelion, estimated at '^lO^^ g (a sphere ^-^150-200 
m across) , still possessed significant cohesive strength. The spine tail — the product of the terminal 
outburst — was a synchronic feature, whose brightest part contained submillimeter-sized dust parti- 
cles, released at velocities not exceeding 30 ms~^. The loss of the nuclear condensation prevented an 
accurate orbital-period determination by traditional techniques. Since the missing nucleus must have 
been located on the synchrone, whose orientation and sunward tip have been measured, we compute 
the astrometric positions of this missing nucleus as the coordinates of the points of intersection of the 
spine tail's axis with the lines of forced orbital-period variation, derived from the orbital solutions 
based on high-quality preperihelion astrometry from the ground. The resulting orbit gives 698 ± 2 
years for the osculating orbital period, which proves that C/2011 W3 is the first major member of 
the expected new, 21st-century cluster of bright Kreutz-system sungrazers, whose existence was pre- 
dicted by these authors in 2007. From the spine tail's evolution, we determine that its measured tip, 
populated by dust particles 1-2 mm in diameter, receded antisunward from the computed position of 
the missing nucleus. The bizarre appearance of the comet's dust tail in images taken only hours after 
perihelion with the coronagraphs on board the SOHO and STEREO spacecraft is readily understood. 
The disconnection of the comet's head from the tail released before perihelion and an apparent activity 
attenuation near perihelion have a common cause — sublimation of all dust at heliocentric distances 
smaller than about 1.8 solar radii. The tail's brightness is strongly affected by forward scattering of 
sunlight by dust. From an initially broad range of particle sizes, the grains that were imaged the 
longest had a radiation-pressure parameter /3 ~ 0.6, diagnostic of submicron-sized silicate grains and 
consistent with the existence of the dust-free zone around the Sun. The role and place of C/2011 W3 
in the hierarchy of the Kreutz system and its genealogy via a 14th century parent suggest that it is 
indirectly related to the celebrated sungrazer X/1106 CI, which, just as the first-generation parent of 
C/2011 W3, split from a common predecessor during the previous return to perihelion. 

Subject headings: comets: general — comets: individual (comet of A.D. 467, X/1106 CI, comet of 
1314, X/1381 VI, C/1843 Dl, C/1880 CI, C/1882 Rl, C/1887 Bl, C/1945 XI, 
C/1963 Rl, C/1965 SI, C/1970 Kl, D/1993 F2, C/2011 W3) — methods: data 
analysis 



1. INTRODUCTION 

The Kreutz system of sungrazing comets (e.g., Kreutz 
1888, 1891, 1901, Marsden 1967, 1989, 2005, Sekanina 
2002a, 2003, Sekanina & Chodas 2004, 2007, 2008) offers 
the ultimate example of an advanced phase of cascading 
fragmentation (Sekanina & Chodas 2007), a process that 
was shown to occur throughout the orbit, including the 
aphelion region (150-200 AU from the Sun). Separation 
velocities acquired by fragments at each fragmentation 
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event generate orbit variations that depend strongly on 
the breakup's heliocentric distance (Sekanina 2002a). At 
a single tidally-triggercd or tidally-assisted event in the 
immediate proximity of perihelion (which takes place in 
the Sun's inner corona), a separation velocity of 1 m s~^ 
causes the two fragments to return to perihelion at times 
that are typically some 80 years apart. The same sep- 
aration velocity at a single nontidal splitting event far 
from the Sun, including the aphelion region, affects the 
orbital period hardly at all, but, depending on the ve- 
locity's direction, can introduce material changes in the 
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other orbital elements: up to ~ ^ Rq (1 Rq = Sun's 
radius or 0.0046548 AU) in the perihelion distance (thus 
allowing some fragments to collide with the Sun's pho- 
tosphere at the next perihelion passage), up to ~5° in 
the longitude of the ascending node and the argument of 
perihelion, and up to ~1° in the orbit inclination. Given 
that the separation velocity can easily reach a few meters 
per second, the resulting effects on the orbital elements 
substantially exceed those by the indirect planetary per- 
turbations considered by Marsden (1967). 

The observed long-term temporal distribution of the 
bright sungrazers does indeed show a tendency toward 
clumping, with the interval between two consecutive clus- 
ters averaging about 80 years (e.g., Marsden 1967, Seka- 
nina & Chodas 2007). In addition, the four major frag- 
ments of comet C/1882 Rl, the brightest known mem- 
ber of the Kreutz system since the beginning of the 18th 
century (Kreutz 1888, 1891) were found to have actual 
orbital periods that increase also with an average step of 
80 years, from ^600 to ^840 years (Sekanina & Chodas 
2007). 

Considering this recurrence cycle, given that the two 
most recent clusters of bright Kreutz sungrazers peaked 
in the 1880s and 1960s, noting that the clusters were pre- 
ceded, two to more than three decades earlier, by precur- 
sor sungrazers, and also recognizing that the arrival rate 
of SOHO Kreutz system comets has been climbing ever 
since the launch of the spacecraft (Sekanina & Chodas 
2007, Knight et al. 2010), the authors predicted in 2007 
that "another cluster of bright sungrazers is expected to 
arrive in the coming decades, the earliest member pos- 
sibly just several years from now" (Sekanina & Chodas 
2007). 

With T. Lovejoy's discovery of C/2011 W3 (Green 
2011), the question has arisen whether this is indeed the 
first major member of the predicted 21st-century cluster. 
The answer depends critically on the accurate determi- 
nation of the comet's orbital period P. If P is about 400 
years or shorter, the comet should be a fragment of one of 
the sungrazers reported in the course of the 17th-century 
or even more recently, and has nothing in common with 
a new cluster. On the other hand, if the orbital period 
is substantially longer than 400 years, then C/2011 W3 
should indeed belong to the new cluster. 

2. PREPERIHELION GROUND-BASED OBSERVATIONS 

When discovered on November 27, C/2011 W3 had 
only ~18 days to reach perihelion (which occurred on 
2011 December 16.0 UT), and there was little hope for an 
accurate determination of the orbital period from preper- 
ihelion data, regardless of their quality. Observing cir- 
cumstances were unfavorable for ground-based imaging, 
since the comet's elongation from the Sun was rapidly de- 
creasing, from 50° at discovery to merely 17°.6 on Decem- 
ber 10, when its position was measured from the ground 
for the last time before perihelion. Still, more than 100 
astrometric positions were obtained during this two- week 
period (Spahr et al. 2011, 2012), the great majority of 
which was sufficiently accurate and mutually consistent 
to be used for deriving a high-quality set of elements, 
except for the orbital period. For example, two sets of 
elliptical elements computed by Williams (2011a, 2011b) 
from 91 observations between November 27 and Decem- 
ber 8 and from 94 observations between November 27 and 



Table 1 

Temporal Coverage of the Head of Comet C/2011 W3 
BY THE SOHO AND STEREO Imaging Instruments^ 





Imaging 


Spatial 


Coverage: 2011 December (UT) 




instru- 


resolu- 






Spacecraft 


ment 


tion'' 


preperihelion 


postperihelion 


SOHO 


C2 


11".4 


15.75-15.96 


16.07-16.22 




C3 


56" 


14.09-15.90 


16.15-18.36 


STEREO-A 


CORl 


7".5 


15.87-15.96 


16.09-16.16 




COR2 


14".7 


15.35-15.92 


16.13-16.66 




HIl 


70" 


12.01-14.92 


16.65-22.37 


STEREO-B 


CORl 


7".5 


15.88-15.96 


16.23-16.45 




COR2 


14".7 


15.35-15.90 


16.38-17.35 




HIl 


70" 


10.81-15.01 


17.87-26.95 


All 






10.81-15.96 


16.07-26.95 



^ Based on the authors' inspection of the images available at these 
websites: http:/ /sohodata.nascom.nasa.sov/cgi-bin/dataj^uery for 
SOHO and http://stereo-sse.nascom.nasa.gov/cgi-bin/iniages. for 
STEREO. 

^ Per pixel or per two binned pixels, as used. 



December 10, gave osculating orbital periods of, respec- 
tively, 376 ± 51 years (leaving a mean residual of ±0".7) 
and 680 ± 64 years (leaving ±0".8). 

3. SPACEBORNE OBSERVATIONS NEAR PERIHELION 
PASSAGE 

Around perihelion, from December 11 to 22, the comet 
was too close to the Sun to obtain any astrometric posi- 
tions from the ground. The comet was, however, exten- 
sively observed with instruments on board several satel- 
lites and deep-space probes. Astrometric data were in 
fact extracted from the comet's images seen in the coro- 
nagraphs and other imaging devices on the SOHO and 
both STEREO spacecraft (Kracht 2011, 2012, Spahr et 
al. 2012). The coverage of the comet's motion by the 
three spacecraft was excellent over this period of time, 
as is apparent from Table 1. We use the standard abbre- 
viations for the relevant instruments: C2 and C3 for the 
two coronagraphs on the SOHO spacecraft and CORl 
for the inner coronagraph, C0R2 for the outer corona- 
graph, and HIl for the first heliospheric imager on either 
of the two STEREO spacecraft. Included in Table 1 is 
the spatial resolution of the instruments, based on the 
information from Brueckner et al. (1995) for SOHO and 
from Howard et al. (2008) for STEREO. The astromet- 
ric data from the measured images of the comet obtained 
by these instruments turned out to be so poor that they 
could not be combined with the ground-based positions, 
often of a subarcsec accuracy, to derive high-quality or- 
bital solutions. The best among the spaceborne data 
are the six positions obtained by Kracht (2011) from the 
C0R2 images of STEREO-B between December 16.49 
and 16.57 UT, shortly after perihelion, but even they are 
not accurate enough to be included in the final iteration 
of the orbit. 

The extensive sets of images taken by the various in- 
struments on board the SOHO and STEREO spacecraft 
prove more useful for examining the comet's dust-tail 
morphology (Sec. 10) and may also be useful for study- 
ing the light curve in the general proximity of perihelion, 
even though many of the CCD frames show the comet's 
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MAJOR CHANGES IN COMET C/2011 W3 (LOVEJOY) ON 2011 DEC. 17-20 




2011/12/17.37 UT 2011/12/18.37 UT 2011/12/19.37 UT 2011/12/20.33 UT 

Figure 1. Some of the earliest postperihelion images of Comet C/2011 W3 taken — at the request of J. Cerny — by J. Ebr, M. Prouza, 
P. Kubanek, and M. Jelmek with the FRAM 30-cm f/10 Schmidt-Cassegrain reflector, a robotic, remotely controled telescope at the Pierre 
Auger Observatory, Malargue, Argentina. Each frame is approximately 11' on a side, corresponding to some 430,000 km in the first image 
and 375,000 km in the last image. North is up, east to the left. The comet was 5°.8 from the Sun during the first exposure, 8°.5 during 
the second, 11°.0 during the third, and 13°.3 during the last one. These images provided the first evidence of the major physical changes 
in the comet's morphological appearance, which culminated with the sudden, complete loss of the nuclear condensation on December 20. 
(Image credit: J. Cerny, Czech Astronomical Society.) 



head saturated. We will not discuss the complex issues 
of the SOHO and STEREO photometry, but would like 
to call attention to a few very preliminary findings from 
our cursory inspection of the images taken with the C2 
and C3 coronagraphs on board SOHO. Most of the in- 
spected images display the saturation artifact known as 
"blooming," which has for similar integration times been 
assumed to measure approximately the amount of excess 
brightness by the number of affected pixels and therefore 
by the length of the overflow streak. Implementation of 
this admittedly oversimplified rule does, however, lead 
to conclusions that arc consistent with the results based 
on firmer, independent evidence addressed later in this 
paper. It appears that the comet's activity became sig- 
nificantly lower very close to perihelion over a period 
of several hours. This apparent attenuation is almost 
certainly a product of intensive sublimation of dust, as 
examined in some detail in Sees. 10 and 11. 

In any case, it seems that the comet's brightness 
reached a maximum about 0.3 day before perihelion, 
followed by rapid fading. Up to three subsequent 
outbursts (some perhaps multiple) may have occurred, 
peaking at, respectively, 0.4, 0.8, and aboutl^ days after 
perihelion, or on December 16.4, 16.8, and ~17.5 UT. 
A rigorous study of the light curve should verify these 
preliminary results. 



4. EARLY POSTPERIHELION OBSERVATIONS FROM 
GROUND, AND SUDDEN TRANSFORMATION OF 
COMET'S APPEARANCE 

Although, astonishingly, the comet was imaged in day- 
light from the ground as early as 1-2 days after perihe- 
lion by several observers (Kronk 2011), the frames con- 
tained no reference stars to derive the comet's astromet- 
ric positions. Among these early postperihelion obser- 
vations was a set of images taken by a group of Czech 
observers between December 17.4 and 20.3 UT, who used 
a robotic, remotely controled telescope, a FRAM 30-cm 



f/10 Schmidt Cassegrain reflector of the Pierre Auger 
Observatory at Malargue, Argentina. While no astrom- 
etry was possible, comparison of the images obtained on 
the four days (Fig. 1) shows major changes in the ap- 
pearance of the comet's head (Cerny 2011). From De- 
cember 17 to 18 the miclear condensation seemed to have 
grown and brightened a little, extending in a broad fan 
in the tailward direction, but otherwise the morphology 
remained essentially the same. Prom December 18 to 
19 the change was more pronounced; even though the 
nuclear condensation remained clearly visible, the quasi- 
parabolic contours of the tail became filled with much 
more material and one can discern a streamer that ex- 
tended for a few arcminutes nearly along the tail axis. 
Should it persist and change its orientation predictably 
with time, such a feature is diagnostic of a sudden, brief 
outburst of dust from the nucleus (Sec. 5). By itself, an 
event of this kind may or may not be part of a cataclysmic 
process that results in the demise of the nucleus. But 
the stunning change in the comet's appearance between 
December 19 and 20 strongly suggests that during this 
episode, portended by the morphological changes during 
the previous days, the nucleus entirely disintegrated. On 
December 20 the nuclear condensation completely disap- 
peared (Fig. 1) and the streamer, usually referred to as 
a spine tail and much longer and more prominent now 
than the previous day, dominated the comet's head and 
near-tail region. 

On December 20 the comet was imaged at Malargue for 
the last time. Although the spaceborne observations with 
the HIl imagers on both the STEREO-A and STEREO- 
B spacecraft continued past this date (Table 1), their 
spatial resolution was not sufficient to show the changes 
detected with clarity in the Malargue images. Fortu- 
nately, the observing conditions on the ground were gen- 
erally improving as the comet's elongation from the Sun 
steadily continued to increase, and after a three-day gap 
earth-based monitoring of the comet's head region re- 
sumed. 
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SPINE TAIL OF COMET C/2011 W3 (LOVEJOY) ON 2011 DEC. 24-2012 JAN. 6 




2011/12/24.74 UT 2011/12/28.73 UT 2012/1/1.72 UT 2012/1/6.74 UT 

Figure 2. The head and adjoining portion of the spine tail of comet C/2011 W3 taken between 2011 December 24 and 2012 January 6 by 
R. H. McNaught with the Uppsala 50-cm f/3.5 Schmidt Telescope, Siding Spring Survey. The fields of view are 48' on a side, corresponding 
to 1.39 million km on the first image and 1.04 million km on the last image. North is up and the east to the left. The Sun is to the upper 
left, but at an angle to the direction of the tail. Clearly seen is the clockwise rotation and a gradual widening of the spine tail, which is 
encompassed by a much broader quasi-parabolic envelope. The surface brightness of the spine tail reaches a flat maximum at some distance 
from the sunward tip; this distance steadily increases with time. The estimated location of the peaJc surface brightness is populated by 
dust particles subjected to solar radiation pressure of /3peak — 0.006-0.007, whose diameter is 400-500 //m at an assumed bulk density of 
0.4 g cm-3. [Image credit: R. H. McNaught, Siding Spring Survey (UA/NASA/ANU).] 



5. FOLLOW-UP POSTPERIHELION GROUND-BASED 
OBSERVATIONS, AND ORIGIN OF SPINE TAIL 

The continuation of the high-rcsohition monitoring of 
the comet was very important for finding out whether 
the loss of the nuclear condensation, so unambiguously 
documented by the December 20 imaging, became indeed 
permanent and irreversible. If the nucleus did indeed dis- 
appear, there would be no definite point to bisect, and 
this condition would seem to thwart any attempt at get- 
ting accurate postperihelion astrometric data and com- 
puting a high-quality set of orbital elements, including a 
well-determined orbital period. 

Starting on December 23, a systematic series of ob- 
servations of what remained of the comet was bcgim by 
McNaught (2012) with the Uppsala 50-cm f/3.5 Schmidt 
telescope of the Siding Spring Survey in Australia. As 
seen from the examples displayed in Fig. 2, the comet had 
nearly the same appearance in all images taken between 
December 23 and January 18, fully confirming the mor- 
phology in the Malargue images of December 20. There 
was no trace of a nuclear condensation, only a gradually 
vanishing "hood" at the stmward end of a faint feature 
with quasi-parabolic contours, on which a much brighter, 
rectilinear, spine-like tail was superimposed. In the ab- 
sence of a better choice, McNaught measured what he 
perceived as possible candidate positions for the tip of 
the tail at its sunward end, arguably the least objection- 
able substitute for the missing nucleus. 

In spite of this handicap, McNaught 's imaging observa- 
tions provide useful information by allowing one to mea- 
sure rather precisely (with a ±1° precision) the system- 
atic clockwise rotation of the spine tail, clearly apparent 
from Fig. 2. Combined with the measurements of the 
streamer on the Malargue exposures from December 19- 
20, the position angle data measured by the authors are 
listed in Table 2. Each of the measured orientations in 
coltunn 2 was compared with a set of calculated posi- 
tion angles of synchronic features, loci of dust particles 
of different sizes subjected to a range of accelerations by 



solar radiation pressure but released from the nucleus at 
the same time (e.g., Finson & Probstein 1968). The dis- 
tribution of the times of particle release, for which the 
synchrones' position angles match exactly the measured 
orientations, have shown a sharp concentration in time, 
with a random scatter of less than ±0.4 day, providing 
evidence that the spine tail was indeed a product of a 
major, fairly brief outburst (or a rapid sequence of out- 
bursts) that peaked on December 17.6±0.2 UT, about 1.6 
days after perihelion; the comet was then 0.144 ±0.012 
AU from the Sun. (For a preliminary report, see Seka- 
nina 2012.) The predicted variations of the spine tail 
until mid-February 2012 are displayed in Fig. 3. 

The event of December 17.6 (be it single or multiple), 
which is possibly identical with the last of the three post- 
perihelion outbursts mentioned at the end of Sec. 3, must 
have begun only a fraction of a day earlier, probably 
about December 17.2 UT, and led to the disappearance 
of the nuclear condensation some 2-2 1 days later and 
to the termination of activity (Sec. 12). Although de- 
tails of this process are unknown at present, it exhibits 
characteristics very similar to those of the cataclysmic 
fragmentation of another sungrazer, comet C/1887 Bl 
(e.g., Kreutz 1901; Marsden 1967, 2005; Sekanina 1984, 
2002a; Sekanina & Chodas 2004). This general conclu- 
sion is stipported by the absence, from December 20 on, 
of any traces of a second tail that would have contained 
freshly ejected dust; such a tail would have preceded the 
spine tail some 5° to 7° in the clockwise direction. 

6. DEVELOPMENT OF A NOVEL TECHNIQUE FOR ORBIT 
DETERMINATION 

There have been attempts to determine the orbital el- 
ements of C/2011 W3 using some of McNaught's mea- 
surements of the tip of the spine tail. The results of 
one such effort have been published by Williams (2011c). 
Although he did not list the residuals, he quoted Mc- 
Naught's remark that "the tip of the 'spine' . . . lies pretty 
much on the front edge of the parabolic hood. It is rather 
less well defined [on Dec. 24] than on Dec. 23. For this 
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Table 2 

Orientation of Spine Tail of Comet C/2011 W3 



Position angle of spine tail 



uaie ^^u 1^ 


measured 


residual^ 


2011 Dec. 19.37 


239° 


+0°15 


20.33 


237 


-0.34 


23.75 


233 


+0.29 


24.74 


231 


-0.48 


26.74 


229 


-0.12 


27.74 


228 


+0.02 


28.73 


227 


+0.11 


29.73 


225 


-0.82 


30.73 


225 


+0.23 


31.73 


224 


+0.27 


2012 Jan. 1.72 


223 


+0.29 


2.73 


222 


+0.37 


3.73 


221 


+0.53 


6.74 


214 


-0.17 


12.45 


52 


-0.34 


16.56 


51 


-0.01 


18.58 


52 


+0.06 



Observer(s) 



Ebr et al.'' 
McNaught 



The difference is: measured minus calculated from a synchronic 
feature generated by dust outburst on 2011 December 17.6 UT. 
b. 



Images reported by J. Cerny. 



reason I cannot be sure just how closely I am measuring 
the same point as on Dec. 23" (Williams 2011c). Even 
though McNaught continued his observations and reduc- 
tions of the positions of the tip well into Jarniary, his 
results from the post-Dec. 24 images have not been pub- 
lished, nor have they been used by Williams to further 
update the orbit. 
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Figure 3. Temporal variations in the position angle (equinox 
J2000.0) of the spine tail, described as a synchronic feature orig- 
inating from an outburst on 2011 December 17.6 UT. The steep 
rate of change between 2012 January 8 and 11 is due to the comet's 
approach to within 1°.5 of the south celestial pole on January 9. 



Fortunately, McNaught communicated the results of 
his continuing observations to one of our colleagues, 
who provided us with the information conveyed (Chesley 
2012). The outcome is 46 astrometric positions of the tip 
of the spine tail, in addition to the three from December 
23-24 already published (Spahr et al. 2012). In the fol- 
lowing we describe a novel technique that we devised to 
exploit this set of McNaught's astrometric observations, 
together with his images of the spine tail. The goal is to 
extract accurate positions of the missing nucleus and to 
employ them in our orbit-determination efforts. 

Every dust particle released from a comet's nucleus 
pursues its own orbit in space, which differs from, and 
is independent of, the comet's subsequent orbit. The or- 
bital deviations are determined by the release time and 
circumstances as well as the magnitude of solar radi- 
ation pressure that the particle has been subjected to 
after release. Radiation pressure forces the particle to 
move in a field of reduced effective gravity compared 
with that acting on the comet. The field intensity is 
described by the ratio /3 between the acceleration 7pi. 
due to radiation pressure and the Sun's gravitational ac- 
celeration 70. If G is the gravitational constant and 
Mq the mass of the Sun, then, at a heliocentric distance 
r, 70 (r) — GMQi/r"^, whereas the acceleration 7pr(7') is 
given as a product of radiation pressure, p(r), and the 
particle's effective cross sectional area for radiation pres- 
sure, ^pr, per its mass, m. Here p(r) = £0/(47rcr^), 
with Cq being the Sun's total radiation energy emitted 
per unit time and c the speed of light, while Ap,- is the 
prochict of the particle's geometrical cross sectional area, 
A, and the efficiency for radiation pressure, Qp^. The ra- 
tio /3 can thus be written in the form (e.g., Sekanina et 
al. 2001) 



Iprjr) _ CQQprA _ 1.15 Qpr 

70 (r) AttcGMqtti px ' 



(1) 



where the last expression on the right applies specifically 
to a spherical grain of a bulk density p (in g cm~^) and 
diameter x (in jim). For large particles (tens of microns 
across and larger) Qpr ~ 1 (the effective cross section for 
radiation pressure is nearly identical to the geometrical 
cross section) and /3 is generally smaller than 0.1. For 
particles released from a comet orbiting the Sun in a 
nearly parabolic path, even this small value of /3 suffices 
to force the dust move in distinctly hyperbolic orbits con- 
cave to the Sun. For strongly absorbing submicron-sized 
grains /? can easily exceed unity (a field with acceler- 
ation by radiation pressure exceeding gravitational ac- 
celeration), and such particles move in hyperbolic orbits 
convex to the Sun. In the special case of /3 = 1, parti- 
cles are subjected to no force and therefore move along 
straight lines. 

As follows from Eq. (1), the ratio /3 is generally in- 
dependent of heliocentric distance. However, the situ- 
ation gets a little more complicated when the comet is 
extremely close to the Sun, as briefly discussed in Sec. 
11. 

It is because of the effects of solar radiation pressure 
that the dust particles of different sizes released during 
a brief outburst line up in the tail along a synchrone, on 
which the nucleus is located at /3 = 0. And because the 
motions of these particles through the tail are indepen- 
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dent of any changes in the comet's motion and behavior 
after the outburst, inchiding any misfortunes that the 
nucleus may then incur, the undisturbed positions of the 
missing nucleus can be recovered from the motion of the 
synchrone, if one can determine the point at which /3 = 0. 
Without an additional constraint the projected position 
of a disintegrated nucleus is increasingly uncertain, but 
primarily in one dimension only, along the synchrone. 

In practice, a synchrone has a finite breadth for a vari- 
ety of reasons: a finite duration of the outburst; being the 
product of a sequence of outbursts; and/or dust particles 
acquiring lateral velocities upon their release. The gen- 
eral tendency for the synchronic features in the dust tails 
of comets is to get broader with time, which gradually in- 
creases the positional uncertainty across the synchrone. 

Over a limited range of fairly low /3 values, the syn- 
chrone is a straight line with very high precision (often 
to better than ±0°.l in the position angle). If the tip of 
the spine tail of C/2011 W3 measured by McNaught is 
located on the synchrone's axis, the coordinates of the 
tip and the spine tail's orientation define the synchrone's 
equation that also fits the undisturbed position of the 
missing nucleus. If [aobs(i), <^obs(i)] are the measured 
equatorial coordinates of the tip of the spine tail and 
p(t) is the tail's position angle at time t, the equatorial 
coordinates of the missing nucleus, [a{t), 5{t)], must sat- 
isfy a condition 
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a{t) ~ ctpbsft) 



cos5obs(t) = tanp(t), 



(2) 



where a{t) and aobs(i) are in hours and 6{t) and (5obs(i) 
in degrees. 

To further constrain the equatorial coordinates a{t) 
and d{t), we recall from Sec. 2 that the quality and con- 
gruence of the astrometric positions from the preperihe- 
lion ground-based observations allowed us to determine 
accurate sets of orbital elements, except for the osculat- 
ing period P. If these early astrometric data are used 
to fit orbits with a number of different forced values 
of the orbital period, the missing nucleus must at any 
given time be located on this line of orbital-period vari- 
ation computed for that time. This condition provides 
a second constraint on the equatorial coordinates [a, S] 
of the missing nucleus, which can readily be found as 
those of the point of intersection of the line of orbital- 
period variation with the line that satisfies the equation 
of synchrone in Eq. (2). This result also provides useful 
information on the orbital period P and determines the 
separation distance of the measured tip of the spine tail 
from the missing nucleus and thereby the critical value 
of /3 for the largest dust grains that are detected at the 
tail's tip. 

Because of the lack of any reliable information on the 
orbital period at the beginning of this exercise, it is 
necessary to search for the solution iteratively. In the 
first approximation, we choose three widely different val- 
ues of the orbital period, P-i, Pq, and Pi, such that 
P_i < Pa < Pi and Pq - P_i = P^ - P^ = APq, derive 
sets of orbital elements from the early astrometry, and 
for each t compute the corresponding topocentric coor- 
dinates [Q!_i(t), [ao{t),6o{t)], and [ai{t),Si{t)]. 
We assume that in the given range of orbital periods 
P{t) both coordinates can be interpolated by fitting, sep- 



arately for each coordinate, a qiiadratic law 

ait) = A{t) + B{t) [P{t)-Po]+C{t) [P{t)-Po]^ 
6{t) = F{t) + G{t) [P{t) - Po] + H{t) [P{t) - Po]^ (3) 

where 



A{t)=aQ{t), 

aijt) - a-i{t) 

= — 2KPa — ' 

a-\{t) + ai{t) 



(4) 



C{t): 



2ao{t) 



2 (APo)2 



and similarly for F{t), G{t), and H{t). After inserting 
for a.{t) and S(t) from Eqs. (3) into Eq. (2) we obtain for 

AP{t) = P{t)-Po: 



[Hit) mat)] [AP(t)]2 + [Git) -Cit)Bit)] APit) 
+Fit)-dobsit) - Cit)[Ait)-aobsit)] = 0, 



where 



<^{t) = 15cos^obs(i) cotp(t). 



(5) 
(6) 



One of the APit) roots of Eq. (5), which must give P(t) 
outside the chosen range of (P_i,Pi), is ignored. The 
other root is used to determine the desired topocentric 
coordinates [ait), Sit)] from Eq. (3). At the end of each 
iteration the values of P(<) are averaged: 



(7) 



where N is the number of images included in the exer- 
cise. It is also verified that scatter of the P values is not 
excessive (not more than a few percent of (P)) and that 
their distribution is essentially random. If (P) differs 
substantially from the starting value of Pq, one should 
contemplate another iteration by replacing this value of 
Pq with (P), tightening the interval APq, and repeating 
the entire procedure, including the determination of the 
time of outburst that controls the position angles of the 
spine tail. The ultimate test is provided by the residu- 
als of the corrected equatorial coordinates of the missing 
nucleus from the new orbital solution. 

7. COMPUTATIONS AND THE FINAL ORBIT 

To keep the number of necessary iterations of the pro- 
posed technique to a minimum, we needed the best possi- 
ble first-approximation set of orbital elements, including 
P. On the one hand, we learnt that the preperihelion ob- 
servations alone were inadequate to provide such a set. 
On the other hand, we found ourselves in an unenvi- 
able situation in regard to the sources of postperihelion 
ground-based data, with no astrometry possible either 
because of the absence of reference stars in images in 
which the comet still possessed a nuclear condensation 
(to bisect in a measuring machine) or because of the loss 
of the condensation (when reference stars were plentiful). 
Under these circumstances, we decided to use some of the 
six positions measured by Kracht (2011) in the STEREO- 
B C0R2 images taken about ^ day after perihelion (Sec. 
3). Combined with the high-quality ground- based posi- 
tions obtained before perihelion, they offered a solution 
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16''58™00^ 16*'57™40^ 16'^7"^0^ 16''57™00^ 
RIGHT ASCENSION 



Figure 4. Schematic illustration of the novel technique developed for determining the corrected topocentric positions of the missing 
nucleus of Comet C/2011 W3 (equinox J2000.0). This example refers to the 15 astrometric positions of the tip of the spine tail measured by 
R. H. McNaught on 2012 January 1 between 17:22:43 and 17:42:25 UT (open circles). Plotted for the first (Point 1) and the last (Point 6) 
positions is the spine tail's orientation (medium-thick lines) as a synchrone that refers to the dust outburst on December 17.6 UT. Projected 
in the direction opposite the direction of the spine tail (broken lines shown again only for Points 1 and 6), each of these lines intersects 
the corresponding line of orbital-period variation, which consists of the predicted topocentric positions of the missing nucleus calculated 
from the preperihelion high-quality observations as a function of the forced orbital period P. Drawn only for the time of Point 1, this line 
of orbital-period variation is shown as a thin line calibrated with a number of P values (in years). The points of intersection, plotted as 
circled dots, are the reconstructed positions of the missing nucleus. We note that the line of orbital-period variation for Point 1 intersects 
the synchrone at a value of P slightly exceeding 700 years. The heavy line is the resulting orbital solution with the ephemeris positions 
at three times, plotted by filled circles. The positions of the spine tail's measured tips are seen to be several arcminutes away from the 
comet's orbit; these distances illustrate the very large magnitude of the introduced errors, if the tip measurements are not corrected. 



acceptable for our purpose, with a maximum residual of 
9", and yielded an osculating period of 785 ± 14 years. 

Based on this finding, we used the same sample of ob- 
servations to compute three sets of orbital elements by 
successively forcing the orbital period to 600, 800, and 
1000 years. At this point we applied the technique de- 
scribed in Sec. 6 to McNaught's 46 astrometric observa- 
tions made between 2011 December 23 and 2012 January 
6, and from Eq. (3) we obtained the first set of predicted 
positions for the missing nucleus. The gist of the proce- 
dure is depicted in Fig. 4 on an example of 15 positional 
measurements obtained by McNaught on 2012 January 
1. The average orbital period resulting from this ex- 
ercise was equal to (P) — 709 ± 6 years. We got ready 
for a second iteration by inspecting the used positional 
data in the orbital run. Having found a solution that is 
much closer to the true orbital period of the comet, we 
were no longer critically dependent on the postperihe- 
lion data points from STEREO-B. They were not used 
in the second iteration, because; the preperihelion astro- 
metric observations alone allowed computer runs with 
new forced values of the orbital period. For the nominal 



run, we chose, in accordance with the developed proce- 
dure, Pq — 709 years, and the resulting new set of orbital 
elements was used to correct the derived position angles 
of the spine tail. We narrowed down the interval of AP 
to ~5fT of (P), or 30 years (i.e., P_i = 679 years and 
Pi = 739 years), and proceeded with the second itera- 
tion of the missing nucleus' positions. Comparing the 
results of the two iteration cycles, we noticed that in 38 
out of the 46 observations the agreement in the derived 
positions of the missing nucleus was 6" or better in ei- 
ther coordinate, and that in the two worst matches the 
differences amounted to 20" and 10", respectively. We 
felt that another iteration was unnc;c;essary because it 
would only lead to changes in the subarcsecond range. 
Consequently, the iterative process was at this point ter- 
minated. 

For the times of McNaught's exposures, the resulting 

topocentric positions of the missing nucleus from the fi- 
nal iteration are listed in columns 2 and 3 of Table 3. In- 
cluded also are the derived positions from the last three 
observations that McNaught made on January 12, 16, 
and 18, although they are too inaccurate to be used in 
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Table 3 



Derived Topocentric Positions of the Missing Nucleus of Comet C/2011 W3 and Related Spine- Tail Data 





Equatorial coordinates 


Distance of 


Acceleration 


Particle 


Residuals'' 


Observation time 






tail's tip from 


parameter, 


diameter, 














t (UT) 


R.A.(2000) 


Decl.(2000) 


nucleus, i^ip 




a^tip (mm) 


R.A. 


Dccl. 


2011 Dec. 23.75224 


h m s 

16 59 21.69 


-39°25'32!8 


1.03 


0.00198 


1.45 


-|-3''5 


-4'8 


24.74267 


16 58 18.93 


-41 44 14.7 


1.40 


0.00202 


1.42 


(+4.4 


-6.1) 


24.74785 


16 58 18.54 


-41 44 57.0 


1.07 


0.00154 


1.87 


+3.2 


-3.8 


26.73669 


16 56 50.13 


-46 40 58.2 


1.20 


0.00107 


2.69 


+0.7 


+1.6 


26.73853 


16 56 50.14 


-46 41 16.6 


0.73 


0.00065 


4.42 


+1.3 


+0.4 


26.74041 


16 56 50.32 


-46 41 38.6 


1.19 


0.00106 


2.71 


+3.8 


-4.0 


27.74339 


16 56 23.95 


-49 20 34.6 


3.04 


0.00221 


1.30 


+0.6 


+1.6 


27.74468 


16 56 24.05 


-49 20 49.1 


2.54 


0.00185 


1.55 


+1.8 


-0.3 


28.73164 


16 56 10.10 


-52 03 48.9 


2.75 


0.00166 


1.73 


(-4.2 


+8.6) 


29.73101 


16 56 09.99 


-54 55 52.1 


2.89 


0.00147 


1.96 


(-5.7 


+9.4) 


30.72663 


16 56 26.64 


-57 54 19.8 


2.79 


0.00122 


2.36 


+3.8 


-2.5 


31.73210 


16 57 00.27 


-61 01 01.5 


4.89 


0.00184 


1.56 


-3.1 


+4.1 


31.73270 


16 56 59.89 


-61 01 05.6 


5.15 


0.00194 


1.48 


(-6.0 


+6.8) 


31.73329 


16 57 00.37 


-61 01 15.4 


5.40 


0.00203 


1.42 


-2.7 


+3.7 


31.73387 


16 57 00.26 


-61 01 21.0 


4.55 


0.00171 


1.68 


-3.7 


+4.8 


31.73447 


16 57 01.36 


-61 01 35.1 


4.70 


0.00177 


1.62 


+4.0 


-2.5 


2012 Jan. 1.72411 


16 57 57.11 


-64 11 34.9 


7.01 


0.00231 


1.24 


(-8.8 


+7.0) 


1.72505 


16 58 00.08 


-64 11 59.9 


7.18 


0.00236 


1.22 


(+10.2 


-6.9) 


1.72599 


16 57 59.62 


-64 12 08.4 


6.99 


0.00230 


1.25 


(+6.7 


-4.3) 


1.72691 


16 57 58.90 


-64 12 15.5 


7.05 


0.00232 


1.24 


+1.5 


-0.6 


1.72786 


16 57 59.43 


-64 12 28.8 


6.52 


0.00215 


1.34 


+4.5 


-2.7 


1.72901 


16 57 58.57 


-64 12 37.8 


6.04 


0.00199 


1.44 


-1.7 


+1.9 


1.72995 


16 57 58.44 


-64 12 47.9 


6.78 


0.00223 


1.29 


-3.0 


+2.8 


1.73089 


16 58 00.61 


-64 13 09.1 


6.65 


0.00219 


1.31 


(+10.7 


-7.3) 


1.73184 


16 58 00.05 


-64 13 17.3 


6.56 


0.00216 


1.33 


(+6.6 


-4.3) 


1.73279 


16 57 59.39 


-64 13 25.0 


6.64 


0.00218 


1.32 


+ 1.8 


-0.8 


1.73395 


16 57 58.98 


-64 13 36.3 


6.83 


0.00225 


1.28 


-1.4 


+1.6 


1.73492 


16 57 59.18 


-64 13 48.4 


6.95 


0.00229 


1.26 


-0.6 


+0.9 


1.73587 


16 58 00.01 


-64 14 03.2 


6.53 


0.00215 


1.34 


+4.4 


-2.7 


1.73682 


16 57 59.46 


-64 14 11.3 


6.19 


0.00204 


1.41 


+0.3 


+0.4 


1.73779 


16 57 59.75 


-64 14 22.9 


7.00 


0.00230 


1.25 


+1.7 


+0.3 


2.72541 


16 59 28.62 


-67 29 52.8 


6.39 


0.00186 


1.55 


-4.9 


+2.7 


2.72634 


16 59 28.49 


-67 30 03.3 


5.40 


0.00157 


1.83 


(-6.2 


+3.5) 


2.72731 


16 59 29.53 


-67 30 18.0 


6.16 


0.00179 


1.61 


-0.9 


+0.6 


2.72826 


16 59 27.73 


-67 30 23.7 


7.11 


0.00207 


1.39 


(-11.8 


+6.4) 


2.72920 


16 59 27.78 


-67 30 34.8 


6.75 


0.00196 


1.47 


(-12.2 


+6.7) 


3.73303 


17 01 45.81 


-70 54 27.0 


7.45 


0.00194 


1.48 


(-20.7 


+6.5) 


3.73398 


17 01 45.84 


-70 54 38.6 


7.76 


0.00202 


1.42 


(-21.3 


+6.7) 


3.73493 


17 01 45.82 


-70 54 50.1 


7.94 


0.00206 


1.40 


(-22.2 


+6.9) 


3.73587 


17 01 44.23 


-70 54 58.8 


7.81 


0.00203 


1.42 


(-30.8 


+9.9) 


3.73682 


17 01 49.66 


-70 55 19.4 


7.59 


0.00197 


1.46 


-5.0 


+1.1 


6.73909 


17 22 41.59 


-81 24 11.4 


8.60 


0.00169 


1.70 


(-28.4 


-5.0) 


6.74001 


17 22 54.11 


-81 24 21.1 


8.65 


0.00170 


1.69 


-2.3 


-3.0 


6.74095 


17 22 39.36 


-81 24 45.4 


8.04 


0.00158 


1.82 


(-37.4 


-15.3) 


6.74189 


17 22 46.40 


-81 24 52.5 


8.51 


0.00168 


1.71 


(-23.7 


-10.4) 


6.74282 


17 22 58.62 


-81 24 55.5 


8.94 


0.00176 


1.63 


+1.6 


-1.6 


12.44813 


4 10 23.94 


-78 38 18.4 


23.29 


0.00340 


0.85 


(-626.6 


-531.2)b 


16.56300 


4 29 30.00 


-65 58 51.9 


12.67 


0.00171 


1.68 


(-98.0 


-125.8) 


18.58108 


4 32 31.89 


-60 35 49.1 


10.75 


0.00143 


2.01 


(-104.9 


-153.3) 



^ The diflcrcncc is: the position derived from measurement minus the position computed from the final orbital solution; the 
residual in R.A. includes the factor cos(Decl.); the positions whose residuals are parenthesized have not been used in the solution. 

^ The starting position of the tail's tip appears to be grossly in error. 



the orbital solution. In column 4 the table lists the an- 
gular distance £tip of the measured tip of the spine tail 
from the nucleus' predicted position. From an ephemeris 
of the December 17.6 synchrone this angular distance 
is converted in column 5 into the radiation-pressure pa- 
rameter /3tip, which in the given range varies practically 
linearly with the angular distance and which, according 
to Eq. (1), is diagnostic of the diameter xtip of the dust 
particles that were detected by McNaught at the tip of 



the spine tail. These are the largest particles in the ob- 
served tail. For an assumed particle bulk density of 0.4 g 
cm~^ (e.g., Richardson et al. 2007), this particle diame- 
ter is given in column 6. Remarkably, Table 3 and Fig. 5 
show that the values of /3tip and Xtip did not vary system- 
atically during the nearly four weeks of McNaught's ob- 
servation, averaging, respectively, 0.00191 ± 0.00042 and 
1.6 ± 0.5 mm. On the other hand, the angular distance 
4ip was increasing nonuniformly, but steadily, with time, 
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resembling the rate of recession of a separated companion 
fragment from the nucleus of its parent comet. 

The tabulated topocentric positions of the missing nu- 
cleus were at this point incorporated into the input file 
for the final orbit determination run. It should be re- 
membered that McNaught's intention was to measure 
positions of the candidates of the tail's tip, not posi- 
tions of extreme points on the tail's axis, which would 
be ideal for our orbit determination efforts. On the days 
on which he took more than one image, it is likely that 
the tail's tip candidates differed from image to image, 
as is abundantly clear from his comments to Williams 
(2011c). And as the spine tail steadily broadened, the 
chance that any measured tip candidate happened to be 
in fact located very close to the axis became increasingly 
remote. One should therefore expect that for only some 
of the derived astrometric positions of the missing nu- 
cleus will the residuals from the final orbit be acceptable 
for the orbit determination run. Also expected is that 
the quality of these positions should deteriorate with the 
spine tail's growing breadth and therefore with time. 



TIME AFTER PERIHELION (days) 
3 13 23 33 
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FROM NUCLEUS OF C/2011 W3 
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J I I L 



J I I L 



• — 
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DATE 2011/2012 (UT) 

Figure 5. Angular distance £tip of tlie spine tail's tip, measured 
by McNaught in his images of C/2011 W3, taken with the Uppsala 
Schmidt Telescope between December 23 and Jajiuary 18, as a 
function of time (top); and the derived diameter of dust particles 
situated at the tip (bottom). The distance increases with time, 
whereas the particle dimensions remain the same. 



The final orbital solution, which includes the rclativis- 
tic effect, is presented in Table 4. It is based on a total 
of 123 observations between 2011 November 27 and 2012 
January 6, all ground-based, which left in either coordi- 



nate a residual not exceeding 5". Of these, 96 arc preper- 
ihelion (from November 27 to December 10) and 27 are 
postperihelion. All postperihelion entries come from Ta- 
ble 3. This means that the position of the missing nucleus 
was successfully recovered from nearly 60 percent of Mc- 
Naught's positions of the tail's tip from images he took 
between Decmcber 23 and January 6. Given the odds 
our approach was facing, we consider this number to be 
quite satisfactory. As expected, the three positions be- 
tween January 12 and 18 turned out to offer unacceptable 
residuals, and the starting position from January 12 was 
clearly incorrect. The residuals of all the positions of the 
missing nucleus are listed in the last two columns of Table 
3, with the rejected positions parenthesized. The resid- 
uals from the preperihelion observations are excellent, 
showing no systematic trends whatsoever, and largely in 
the subarcsecond range. The two most significant contri- 
butions, with images in both cases obtained robotically, 
were from the Pierre Auger Obsca'vatory, Malargue, Ar- 
gentina (code 147, 42 accepted positions) and from the 
Remote Astronomical Society Observatory, Mayhill, New 
Mexico (code E03, 21 accepted observations); combined 
they account for nearly two-thirds of all used preperihe- 
lion positional data. 

Out of a total of 46 positions available to us from as- 
trometric measurements of the images taken on board 
the SOHO and STEREO spacecraft, none satisfied our 
adopted cutoff of ±5" for the residuals. As expected 
from our discussion in Sec. 3, the least discordant were 
again the measurements by Kracht (2011) from the six 
C0R2-B images of Dec. 16, which left residuals of up to 
27" in right ascension and up to 15" in declination. The 
24 positions from the STEREO-B images between Dec. 
10 and 14 showed surprisingly small residuals in decli- 
nation, of up to 12", but strongly systematic, entirely 
unacceptable residuals of up to 60" in right ascension. 
The 8 data points from the SOHO images on Dec. 14 
had the opposite problem; 5 of them had residuals less 
than 10" in right ascension, but all were off by 40" to 
87" in declination. The least satisfactory was a set of 8 
positions from STEREO-A (hsted by Spahr ct al. 2012), 
which had residuals scattered wildly from —94" to +97" 
in right ascension and from —78" to +4" in declination. 

From the orbital period in Table 4 we conclude that 
comet C/2011 W3 is the first major member of the new, 
21st-century cluster of bright Kreutz sungrazers whose 
existence was predicted by the authors of this paper in 
2007 (Sec. 1). The comet cannot be the return of a frag- 
ment of any of the sungrazing comets observed since the 
17th century, contrary to speculations based on prelim- 
inary sets of orbital elements. The perihelion distance 
agrees closely with the perihelion distances of the bright 
sungrazers C/1843 Dl, C/1880 CI, and to a lesser degree 
C/1963 Rl, but all three angular elements have values 
unlike any of the known bright members of the system.^ 

^ As a matter of historic curiosity, it should be pointed out that 
the first orbit for the sungrazer C/1945 XI, calculated by Cun- 
ningham (1946a, 1946b). was in fact somewhat similar to the orbit 
of C/2011 W3 {uj = 50°.9, n = 322°.3, i = 137°.0, q = 0.0063 
AU), but it was derived from only crudely determined positions 
of the comet. The photographic plates with the comet's images 
were properly measured and reduced only in 1952 (Marsden 1967) 
and none of Marsden's (1989) resulting orbits shows any obvious 
resemblance to that of C/2011 W3. 
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Table 4 

Final Set of Orbital Elements Adopted for Comet C/2011 W3 and 
Integration One Revolution Back in Time (Eq. J2000.0) 



Orbital clement 



Current apparition 



Previous return'^ 



Epoch of osculation (ET) 
Time of perihelion passage, tn (ET) 
Argument of perihelion, u) 
Longitude of ascending node, Q 
Orbital inclination, i 
Perihelion distance, q (AU) 
Orbital eccentricity, e 
Osculating orbital period, P (yr) 



2011 Dec. 25.0 
2011 Dec. 16.011810 ± 0.000040 
53°.5103 ± 0°.0020 
326°.3694 ± 0°.0027 
134°.3559 ± 0°.0012 
0.00555381 ± 0.00000007 
0.99992942 ± 0.00000014 
698 ±2 



1329 Jan. 6.0 
1329 Jan. 4.9336 
52°.815 
325°. 276 
133°.811 
0.0059198 
0.9999237 
684 



^ Dates in the Julian calendar. 



This suggests that there may exist yet another subcate- 
gory of bright sungrazers that has never been considered 
in the evolutionary models of the Kreutz system. Table 
4 also shows the results of integrating the orbit back one 
revolution about the Sun, to the early 14th century, sug- 
gesting the presence of fairly minor perturbations. The 
related issues and implications of the comet's orbit are 
addressed in Sec. 13. 

8. PARTICLE VELOCITIES IN SPINE TAIL AND ITS 
ANOMALOUS BRIGHTNESS PROFILE 

Besides the clearly apparent clockwise rotation of the 
spine tail, both the Malargue images in Fig. 1 and the 
Siding Spring images in Fig. 2 show a quasi-parabolic 
envelope, in which the spine tail is immersed and which, 
based on its approximate axial orientation, may have 
originated from the previous outburst that peaked on 
about December 16.8 UT or 0.8 day after perihelion (Sec. 
3). While the spine tail dominates the images brightness- 
wise, the previous outburst was apparently more violent, 
as suggested by a much greater breadth of the envelope 
and, consequently, much higher ejection velocities of the 
released dust. Indeed, if the width of the spine tail and 
the envelope are both due entirely to lateral velocities of 
the particles at the time of release, these velocities come 
out to be typically up to 20 30 m s^^ for the spine tail, 
but up to at least 150 m s^^ for the envelope. Whereas 
the latter velocities are fairly typical of microscopic dust 
ejected due to an interaction with outflowing gas, the 
former are so unusually low that they suggest that the 
entire residual mass of the nucleus simply collapsed. 

One may question whether the episode that led to the 
formation of the spine tail can at all be called an out- 
burst. It may be more appropriate to refer to it as a 
cataclysmic fragmentation event or terminal collapse or 
complete disintegration, but we continue to use all these 
terms interchangeably even after we address the physical 
nature of the process in Sec. 12. 

The distribution of surface brightness along the spine 
tail, which displays a broad maximum at some dis- 
tance from the sunward tip in all McNaught's images, 
is positively peculiar. The gradual rise in the near- 
head brightness with increasing distance from the nu- 
cleus is anomalous in that it contradicts a typical power 
law for the particle-size distribution of dust comets. 
In a low-ejection-velocity approximation, the Finson- 
Probstein (1968) approach provides the following expres- 
sion for the variation of the surface brightness I as a 



function of the radiation-pressure acceleration parame- 
ter 13 along a synchronic feature observed at a time fobs 
and made up of dust particles released from the nucleus 
at a time tobs - 



r: 



pcpia) (3 



k-4 



(8) 



where p is the geometric albedo of the particles, (l){a) 
their phase function at a phase angle a [normalized to 
(j){0°) = 1], x~^dx is the differential distribution function 
of particle diameters x, r is the heliocentric distance, and 
9 is a Jacobian that converts the radial and transverse 
rectangular coordinates, ^, rj, into /3 and r 







dp 


dr 


dr] 


dr] 


dp 





(9) 



At a geocentric distance A, the coordinates ^ and r] are 
projected onto the tangent plane of the sky at the comet's 
nucleus and reckoned from it, ^ pointing away from the 
Sun and r] being normal to it in either direction, as only 
the absolute value of 9 is used in Eq. (8). 

We have already mentioned that at small values of (3 
(say, f3 < 0.1), the population of particles released at a 
time fobs — T and accelerated at a rate p is located at 
an angular distance i from the nucleus, which is propor- 
tional to P, 

m = Kp, (10) 

where K is a constant. Similarly, for small values of P, 

d£,/dp (X P-A and dri/dp oc P-A, while d£,/dT(xA and 
drj/dr (X A, so that for such large particles 



9(^,r,;/3,r)(x/3-A2 



and 



DC 



r2A2 



(11) 



(12) 



Since the exponent k for the size-distribution function of 
cometary dust is typically in the range of 3.5 to 4.0 (e.g., 
Sitko et al. 2011), the surface brightness along the spine 
tail of C/2011 W3 is expected to drop with increasing 
distance £ from the nucleus with some power oi £, usually 
between and From Eq. (12) it is obvious that 
the condition d{logI)/d{logi) in the part of the 
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Figure 6. Preliminary light curve of comet C/2011 W3 based on visual and CCD total brightness estimates made from the ground. 

The magnitudes plotted have been normalized to 1 AU from earth with the inverse square power law and to a zero phase angle with the 
"compound" Henyey-Greenstein law, as modified by Marcus. Personal and instrumental effects have been corrected to the limited degree 
possible. The comet is shown to have been much brighter after perihelion than before it. 



spine tail between the tip and the peak surface brightness 
implies a condition 

fc>5. (13) 

Such an extremely steep size distribution function can 
only be explained by cnmibling of the largest particles 
(say, X > 1 cm) into progressively smaller ones. Crude 
visual inspection of McNaught's spine-tail images sug- 
gests that the surface brightness attains the peak near 
/Spoak — 0.006-0.007, or at particle diameters of ^400 to 
500 yum, independent of the observation time. It is there- 
fore possible that the fragmentation process essentially 
terminated once fragments of large debris were reduced 
to sizes of -^100 /xm. 

9. LIGHT CURVE, TAIL BRIGHTNESS, AND TOTAL MASS 
ESTIMATE FOR RELEASED DUST 

It is hoped that eventually a detailed light curve will 
be available for comet C/2011 W3. At this time, how- 
ever, only fragmentary information has been published. 
With the brightness data collected primarily from Spahr 
et al. (2011, 2012) and Green (2012a, 2012b), Figure 6 
presents a preliminary light curve, a plot against the he- 
liocentric distance of the magnitude, H'^, normalized to 
1 AU from Earth with the inverse square power law of 
the geocentric distance and to a zero phase angle with 
the "compound" Henyey-Greenstein law, as modified for 
cometary dust by Marcus (2007). The plotted magni- 
tudes were all identified by the observers as total, either 
visual or CCD, and were corrected by us for personal and 
instrumental effects to the limited degree possible and 
referred to the visual photometric system of an average 



unaided eye. The phase-effect correction is important es- 
pecially for the postperihelion observations, for which the 
phase angle was varying between 90° and 130° and the 
effects of forward scattering of sunlight were significant. 

The comet was much brighter after perihelion than be- 
fore, unquestionably due to its cataclysmic fragmenta- 
tion (Sees. 5 and 8). Figure 6 does not explicitly show 
the outbursts that are mentioned in Sec. 3, but this is 
hardly significant given the poor coverage of the light 
curve near the Sun and low accuracy of the data. From 
Fig. 6 it appears that to a first approximation, one can 
fit both branches of the light curve with the traditional 
power law r~", with n equal to 6.9 ± 0.4 before perihe- 
lion and 5.5 ± 0.2 after perihelion, and with the intrinsic 
magnitude (at 1 AU from both the Sun and earth), Hq, 
equal 15.0 ± 0.3 before perihehon and 10.0 ± 0.2 after 
perihelion. It is possible that without the outbursts the 
postperihelion light curve would have gotten fiattened at 
heliocentric distances below ^0.2 AU. 

The primary purpose for incorporating the light curve 
into the scope of this paper has been our intention to use 
it for estimating the total mass of dust involved in the 
fragmentation process culminating on December 17.6 UT 
and therefore the residual mass of the nucleus at breakup. 
If the visual brightness at the time were due entirely to 
scattering of sunlight by dust released from the disinte- 
grating nucleus into the atmosphere, one could readily 
compute the total cross sectional area of the particulate 
material, once a value for the geometric albedo of the op- 
tically thin cloud is adopted. In practice, unfortunately, 
sodium atoms radiating in the doublet near 5900 A are a 
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Table 5 

Upper Limit on Total Amount of Dust prom Fragmentation of Comet C/2011 W3 



Upper limit on 
particle size'', 
a^max (cm) 


Toial mass 


of (lusi, (.V-lj 




Ellrctivc- . 


liaiiKU.ei'. 


Pma.v (km) 


k = 3.5 


k = 3.8 


fe = 4.1 


k = 3.5 


k = 3.8 


i fe = 4.1 


10 


0.9 X W^* 


5.5 X 10^2 


7.7 X 10^1 


0.76 


0.30 


0.15 




3.0 X 10" 


9.0 X 10^2 


8.2 X 10^1 


1.12 


0.35 


0.16 


10^ 


9.3 X 10" 


14.5 X 10^2 


8.6 X 10" 


1.65 


0.41 


0.16 



^ Lower limit to particle diameters is always a^min = 0.1 jum, the bulk density pdust = 0.4 g cm 



strong contributor to the visual brightness of SOHO sun- 
grazing comets (Bicscckcr ct al. 2002), so that the hght 
curve provides us with only an upper limit to the cross 
sectional area of the dust particles and indirectly to the 
mass of the disintegrated nucleus. 

An independent approach to determining the amount 
of released dust is to estimate the total visual brightness 
of the comet's tail at a time sufficiently long after the 
fragmentation process has been completed. Because the 
tail includes all dust released by the comet since its per- 
ihelion passage, the result in this case is an estimate of 
the mass of the nucleus at, or very shortly after, perihe- 
lion. This result would represent a lower limit: not only 
would all sodium have sublimated away and ionized long 
before such an observation, but by this time the surface 
brightness of some of the released debris would have al- 
ready dropped below the detection threshold of visual 
observers. 

The immediate goal of either of the two approaches is 
the determination of Heft, the effective visual magnitude 
of the estimated light scattered by the debris, normal- 
ized to the heliocentric and geocentric distances of 1 AU 
with an inverse square power law, to a zero phase angle 
with the Marcus (2007) formula, and to the photometric 
system of an average unaided eye. Expressed in km^, 
the effective cross sectional area of the dust in the cloud, 
^dust) is given by an equation 



Xaust = — X lO^-o-^^^e, 
P 



(14) 



where p is the geometric albedo of the dust, for which 
we use a value of 0.04 (e.g., Lamy et al. 2009). For the 
Sun's visual magnitude we adopt in Eq. (14) a value of 
—26.65 (corresponding to the V magnitude of —26.75). 
There is no explicit phase function involved because Hcs 
has already been corrected for this effect. 

To be able to estimate the total mass of dust that corre- 
sponds to the derived cross sectional area, one has to as- 
sume the following: (i) an upper limit on diameters of the 
particles, Xmax, (ii) the lower limit, Xmin, (iii) their size 
distribution function between x^i^ and Xmax, which, as 
in the preceding section, is given by a power law x~^dx, 
and (iv) their average bulk density, pdust- The relation 
between the total mass of dust, A^dust, and the cross 
sectional area, Xjust) is given by 



M 



dust 



2(fc-3) 



3(4 - k) 
when 3 < fc < 4, and by 

2(fc- 



Pdust ^dust X'^i^X 



k-3 4-k 
min ■''max 



^dust = 



3) 



3(A;-4) 



©Pdust -'^dust 



(15) 



(16) 



when A; > 4. The coefficient 9, usually near unity, is 
given by an expression 



e 



('^min/'^max) 



\k-i\ 



1 (^min/ ^max)'^ ^ 



(17) 



We use Pdust = 0.4 g cm^"^ (e.g., Richardson et al. 2007), 
which is probably a good estimate for larger grains, but 
must be an underestimate for microscopic dust. As for 
a^min, particles smaller than 0.1 jim in diameter were de- 
tected in comets, e.g., in the atmosphere of IP/Halley 
with the mass spectrometers on board the Giotto and 
VEGA spacecraft, but their contribution to the total 
mass of dust was found to be low (Utterback & Kissel 
1990, 1995; Sagdeev et al. 1989), so that adopting a lower 
limit of Xmin = 0.1 jim seems to be justified. 

The exponent k of the power law and the upper limit 
a^max were in this exercise varied. In conformity with 
the findings for a large number of comets, k has been 
constrained to a range of 3.5 to 4.1 (e.g., Sitko et al. 
2011), while a;inax has been selected to range from 10 cm 
to 10 m (e.g., Harmon et al. 2011). 

We can now compare the approach based on the light 
curve (providing an upper limit) with that based on 
the tail's brightness (offering a lower limit). Turn- 
ing to Fig. 6, we find that at the time of cataclysmic 
fragmentation on December 17.6 UT, at a heliocentric 
distance of rrci = 0.144 AU, the corrected magnitude 
H'^(rj:ci) — —1.6, which implies that (i/off)max = +2.6 
and an upper limit (Xdust)max = 3.5 x 10® km^. The 
estimated uncertainty is at least ±30 percent. For nine 
different combinations of the parameters k and Xmax; the 
resulting upper limits on the mass (Aldust)max and the 
effective diameter Pmax are listed in Table 5. 

The source of the best information that we have been 
able to find for the brightness of the dust tail of C/2011 
W3 at this early phase of investigation is an attempt 
undertaken by Seargent (2011) on 2011 December 24.69 
UT. Using an out-of-focus method of brightness estima- 
tion for comparison stars, he determined the total bright- 
ness of the dust tail in the field of view (of an estimated 
diameter of ^3°) of his 25x100 binocular telescope at 
three angular distances £ from the head that he esti- 
mated at magnitude 4.8. The results of Seargent's ob- 
servation are presented in columns 1 to 3 of Table 6. 
The four subsequent columns provide information on the 
particles and the expected tail orientation at each dis- 
tance from the head for two assumed times of release. 
The geocentric and heliocentric distances, the phase an- 
gle, and the resulting effective visual magnitude Hes{£), 
all of which are — to the given precision — indepen- 
dent of the choice for the time of release, are listed in 
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Total Brightness op Dust Tail of Comet C/2011 W3 on 2011 December 24.69 UT (Estimates by D. Seargent) 



Angular 
distance i 
from head 


Visual 
magnitude 
estimate 


Observed 
tail 
width 


Parameter 


^ for (t,AY 


Position angle for (trel)'' 


Distance (AU) from 


Phase 
angle 


Effective 
magnitude 

i?aff(^) 


+0.4 day 


+1.6 days 


+0.4 day 


+1.6 days 


Earth 


Sun 


4° 


4 


0°.5 


0.08 


0.36 


237°.9 


231°.7 


0.64 


0.51 


117° 


6.8 


8 


4.5 


1 


0.17 


0.71 


238.4 


231.9 


0.62 


0.56 


113 


7.0 


23 


5 


2 


0.50 


2.09 


239.8 


232.4 


0.58 


0.72 


97 


6.6 


28'' 


c 




0.62 


2.58 


240.1 


232.6 


0.58 


0.77 


93 





^ Assumed time of dust release reckoned from the time of the comet's perihelion passage. 
^ Length of the tail seen by Seargent with the naked eye. 

^ According io Scargent's report, over its last 5° the tail was very faint and difficult with the naked eye and could not be traced in 
the 25x100 binocular telescope. 



columns 8 to 11. The fundamental conclusion from the 
tabulated values of Hcff{i) is that the total brightness in 
each band of the tail was virtually independent of the 
distance from the head and equal to 6.8 ± 0.2: the sur- 
face brightness was decreasing with distance, but the tail 
width was increasing. And since Seargent could detect 
no tail with his binocular telescope beyond ^23° from 
the head, we assume in our conservative estimate that it 
ended there. The whole tail was then brighter by a fac- 
tor of or 2.2 magnitudes, than the average in each 3° 
band, which leads to a lower limit of (-ffcff)min = 4.6, or 
2 magnitudes fainter than indicated by (-ffeiiOmax from 
the light curve. Consequently, the lower limit on the 
total cross sectional area of the comet's dust released 
after perihelion is (X(just)min = 0.6 x 10® km^. To ob- 
tain lower limits on the mass and diameter estimates, 
the values in Table 5 should be multiplied by 0.16 and 
0.54, respectively. 

It appears that the derived lower limit is relatively se- 
cure. Seargent's value for the brightness of the head is 
about 0.9 magnitude below the light curve in Fig. 6 (sug- 
gesting a personal magnitude scale difference) and our 
truncation of the tail's length also adds a few tenths of a 
magnitude, because it implies a loss of light in the binoc- 
ular magnitudes relative; to those obtained with the un- 
aided eye. A fairly narrow range resulting from the two 
entirely independent approaches applied is encouraging 
and may indicate that upon the terminal breakup all 
atomic sodium escaped very rapidly, leaving hardly any 
lasting signature in the comet's light curve, given that at 
0.144 AU from the Sun the Na photoionization lifetime 
is less than 1 hour (Huebner, Keady, & Lyon 1992; also: 
Cremonese et al. 1997, Combi, DiSanti, & Fink 1997). 

Because of the unacceptably large mass estimates, 
greatly exceeding lO'^'^ grams, we suggest that the ex- 
ponent k of the size-distribution power law in Table 5 
could not possibly be near or bcilow 3.5. Howcvcir. there 
are fairly good odds that the nucleus could have been 
as large as 150-200 meters in diameter and that it could 
have had a residual mass as much as -^10^^ grams shortly 
after passing through perihelion. 

10. PROPERTIES OF DUST IN THE TAIL BEFORE AND 
AFTER PERIHELION 

Examining eleven of the brighter SOHO Kreutz-system 
comets with prominent tails, all on their way to perihe- 
lion between 1996 and 1998, Sekanina (2000) noticed that 
the narrow and nearly straight tails deviated strikingly 
from the antisolar direction. The highest solar radiation- 



pressure accelerations, to which the dust in the tails was 
subjected, was found to have been (3 ^ 0.6, suggesting 
the presence of a population of submicron-sized grains di- 
electric in nature, most probably silicates (e.g., Sekanina 
et al. 2001 and the references listed there in Sec. II. A. 2). 
The sampled comets also showed rather consistently that 
the peak production of dust occurred some 20 -30 Rq 
from the Sun, about 1 day before perihelion, and that 
activity essentially terminated shortly afterwards. De- 
tailed modeling of one of these objects showed that its 
tail was a syndyname of /3 = 0.6 far from the head, but 
a synchrone referring to a release time of 0.8 day be- 
fore perihelion (19 Rq from the Sim), near the head. A 
sharp bend or knee, but no gap, separated the two parts 
of the tail. The author suggested that the time of the 
abrupt termination was possibly nuclear-size dependent 
(Sekanina 2000). 

To learn more about this possible relation and about 
the dust-emission pattern of comet C/2011 W3, we have 
examined the properties of its tail both before and after 
perihelion. Two images taken by the C3 coronagraph 
on board the SOHO spacecraft shortly before perihelion, 
on December 15.504 and 15.796 UT, are reproduced in 
Fig. 7. Cursory inspection reveals a tail whose curva- 
ture is slight but increases with time, with no obvious 
knee. Two syndynames are plotted for each image in the 
right-hand side panel: the one on the left corresponds to 
/3 = 2.5, the other to /3 = 0.6. The two radiation-pressure 
accelerations have been chosen because they are of par- 
ticular significance for cometary dust (e.g., Sekanina et 
al. 2001): the first is the highest acceleration known to 
affect dust in comets and is typical of strongly absorbing 
grains (such as carbon-rich, organic material), whereas 
the second, as mentioned above, is characteristic of the 
peak acceleration on dielectric grains (such as silicates). 
In both cases the particles involved are in the submicron- 
size range. Superposition of the syndynames on the im- 
ages in Fig. 7 shows for both observation times that the 
tail lies largely between the two syndynames and that it 
is slightly less curved than either of them. A conclusion 
from these images alone is that C/2011 W3 was releasing 
both dielectric and absorbing dust on its way to perihe- 
lion. The main difference between the two syndynames 
in Fig. 7 is that at a given distance from the nucleus, the 
dust on the syndyname of /3 = 0.6 was released much 
earlier than that on the syndyname of /3 = 2.5. 

It is instructive to compare the comet's appearance in 
the two preperihelion images with very different views 
offered shortly after perihelion by the SOHO C2 corona- 
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SOHO IMAGES AND MODELS OF DUST TAIL OF C/2011 W3 ON DEC. 15 
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Figure 7. Comet C/2011 W3 with its dust tail siiortly before reacliing periiielion, as imaged by the C3 coronagraph on board the SOHO 
spacecraft. The upper image was taken on December 15.504 UT, or 0.508 day before perihelion, the lower image on December 15.796 UT, 
or 0.216 day before perihelion. In the panels next to the images the observations are compared with two syndynames; the one for /3 = 2.5 is 
slightly to the left of the other, which refers to /3 = 0.6. The syndynames are calibrated by the times of dust release, reckoned in days from 
the time of the comet's perihelion passage (the negative numbers indicate days before the comet's perihelion). The panels also provide the 
scale and orientation of the images. We note that in both images the tail displays a slightly lesser curvature than the syndynames, but 
lies generally between them, which implies that the comet released submicron-sized dust of both dielectric and strongly absorbing nature. 
(Image credit: ESA/NASA/LASCO consortium.) 



graph in combination with the CORl and C0R2 corona- 
graphs on board the STEREO-A and STEREO-B space- 
craft. The three selected images are of course merely 
snapshots of continuous changes in the comet's figure 
that are fully revealed only by the dynamic, time-lapse 
imaging that involves the entire data set. Nevertheless, 
given the spatial distribution of the three spacecraft, the 
stereoscopic quality of the gained information makes even 
snapshot views extremely valuable. 
Figure 8 shows the comet's head and a part of its tail 



in a frame taken with the C2 coronagraph on Decem- 
ber 16.117 UT, or 0.105 day after perihelion. The bright 
tail, to the south-southeast of the Sun, is entirely dis- 
connected from the head. The panel to the right of the 
image, showing the complete syndyname /? — 0.6, sug- 
gests that, contrary to the observation, the tail should 
have reappeared to the northwest of the occulting disk. 
The comet's head, although clearly saturated and dis- 
playing some "blooming" , is essentially stellar in appear- 
ance. This means that no detectable postperihelion tail 
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Figure 8. Appearance of eomet C/2011 W3 and its dust tail in an image taken with the C2 coronagraph on board the SOHO spacecraft on 
December 16.117 UT, or 0.105 day after perihelion. The tail, to the south-southeast of the Sun, is seen to be completely disconnected from 
the comet's head, to the west of the Sun. The head's image is saturated (with some blooming being apparent) but with no tail extension. 
The tail lies entirely between the syndynames /9 = 0.6 and 0.8, as shown in the panel to the right. The first of the two syndynames is 
plotted thicker and is calibrated by the times of dust release, reckoned in days from the time of the comet's perihelion passage, pinpointing 
the locations of particles released at perihelion and 0.1, 1, and 2 days before perihelion. The dot on the segment of the syndyname /3 = 0.8 
identifies the location of dust released 1 day before perihelion. We note that the entire northwestern branch of the tail is missing. The tail's 
brightness is largely determined by forward scattering of sunlight by dust, as the phase angle increases toward the bottom of the image. 
The scale and orientation of the image are the same as those of the panel. (Image credit: ESA/NASA/LASCO consortium.) 



developed by the time the image was taken, about 0.1 
day after perihelion. 

An example of the peculiar appearance of C/2011 W3 
is presented in Fig. 9, which shows the comet imaged 
with the C0R2 coronagraph on board the STEREO-A 
spacecraft on December 16.246 UT, or 0.234 day after 
perihelion. Unlike in the SOHO's C2 coronagraph, the 
comet now consists of three components: the head — 
with a short wisp of a new tail that must have begun 
to develop before 0.2 day after perihelion — to the east- 
northeast of the Sun and two separate branches of the old 
tail, one to the southeast and the second to the northwest 
of the Sun. The southeastern branch is relatively sharp 
but faint, a far cry from its luster in Fig. 8. It can be 
matched with the syndyname j3 = 0.6, but the resolution 
is poor as the syndynames along this tail are "crowded" . 
The northwestern branch is blob shaped. Its first view 
in this C0R2 coronagraph coincided approximately with 
the re-appearance of the comet's head from behind the 
occulting disk, nearly 3 hours after perihelion. In pro- 
jection onto the plane of the sky, this branch may un- 
fortunately have been superposed on top of a weak but 
broad coronal mass ejection. Because of this interfer- 
ence, the tail's exact contours are hard to establish, but 
the syndyname P = 0.6 seems to be again involved with 
the feature. The blob terminates just before crossing the 
synchrone for a release time of 0.1 day before perihelion. 
Again, there is no tail in the areas corresponding to re- 



lease times near perihelion. Neither of the two branches 
becomes obvious in any postperihelion image taken with 
the CORl coronagraph of the STEREO-A spacecraft. 

In a long series of images taken with the CORl and 
C0R2 coronagraphs on board the STEREO-B spacecraft 
during much of the first day after perihelion, the look of 
the comet with its tail is downright bizarre. In fact, 
hours before the comet's head emerged from behind the 
occulting disk to the west-southwest, a second branch of 
the tail began to show up as a steadily growing sharp 
spike to the northeast of the Sun, joining the southwest- 
ern branch that had thrived since preperihelion times. In 
addition, the comet's head, after it emerged, was, just as 
in Fig. 9, disconnected from either of the two branches 
of the tail. Overall, therefore, the comet again consisted 
of three discrete components, as seen in Fig. 10. The 
displayed image was taken with the C0R2 coronagraph 
on December 16.517 UT, or 0.505 day after perihelion. 
The syndyname /3 — 0.6 provides an excellent fit to both 
branches of the tail simultaneously, although the syndy- 
names are again fairly "crowded" along much of the tail. 
This image shows that the northeastern branch extends 
to a point that is populated by dust released ^0.1 day 
before perihelion, in agreement with the result from Fig. 
9. Careful inspection of the comet's head reveals its elon- 
gation similar to that detected in the STEREO-A image. 

To extract information on the dust particles that pop- 
ulate the tail, we have examined more closely the or- 
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Figure 9. Peculiar appearance of comet C/2011 W3 and its dust tail in an image taken with the COR2 coronagraph on board the 
STEREO-A spacecraft on December 16.246 UT, or 0.234 day after perihelion. The tail has two branches located, respectively, to the 
southeast and northwest of the Sun, while the comet head's site is nearly to the east of the Sun. The tail's southeastern branch is faint 
because of backscatter of sunlight, while the blob-shaped northwestern branch appears to be at least in part contaminated by a relatively 
weak, but broad coronal mass ejection. The panel to the right identifies the locations on the syndyname /3 = 0.6, which are populated 
by dust released from the nucleus at six different times between 8 days before perihelion and 0.03 day after perihelion. Three additional 
release times — 0.1, 0.3, and 0.7 day before perihelion — are identified by the synchrones drawn by the broken curves. The panel also 
shows the area of the tail's blob to the northwest of the Sun and short segments of the syndynames for /3 equaling 0.5 and 0.8. No tail is 
detected for release times between 0.1 day before perihelion and about 0.1 day after perihelion, but the head's image is elongated in the 
direction away from the Sun, suggesting that a postperihelion tail already began to develop during a couple of hours before the frame was 
taken. (Image credit: NASA/SECCHI consortium.) 



bital properties of dust particles along the syndyname 
(3 = 0.6 with the help of Table 7, which lists their peri- 
helion distance, the eccentricity, and the perihelion time 
as a function of /3 and the time of release from the comet. 
The numbers along this syndyname in Figs. 8-10 are the 
times of release in days from the comet's perihelion time. 
While it is rather clear that the bright branch of the tail 
to the south-southeast in Fig. 8 is identical with the fairly 
pale southeastern branch in Fig. 9 and with the gleam- 
ing southwestern branch in Fig. 10, the reason for the 
enormous differences in brightness is not obvious. The 
orbital elements of the dust particles in Table 7 allow 
one to calculate the geometry for the times of the three 
images and to find out that the culprit is forward scatter- 
ing of sunlight. For example, the particles in the bright 
portion of the tail in Fig. 8, located on the syndyname 
/3 = 0.6 and released 2 days before perihelion, were at 
the time of observation nearly 7 Rq from the Sun and 
their phase angle was 113°, while for the particles that 
left the nucleus 1 day before perihelion, although they 
were only 3.6 Rq from the Sun, the phase angle was 90° 
and for those released 0.7 day before perihelion and lo- 
cated near the edge of the occulting disk, the phase angle 
was only 70°, the reason why this part of the tail looks 
faint. For the same reason, the brightness of the tail's 
southeastern branch in Fig. 9 cannot compare with that 



of the southwestern branch in Fig. 10 even though it is 
the same feature: at the location of particles on the syn- 
dyname /? = 0.6 released 5 days before perihelion, in the 
middle of this tail, the phase angle is 55° in Fig. 9, but 
109° in Fig. 10. 

It is further noted from Table 7 that particles on that 
same syndyname that were released earlier than 2.5 days 
before perihelion did not reach their perihelion points un- 
til after December 16.5 UT and, at the times the images 
in Figs. 8-10 were taken, were still on their way toward 
the Sun. By contrast, particles released later than 0.5 
day before perihelion reached their perihelion points no 
later than about December 16.1 UT and were therefore 
in all three images already moving away from the Sun. 
This is important because the relevant perihelion dis- 
tances are very small. Table 7 shows that on the syndy- 
name l3 = 0.6 this minimum distance drops from 2.37 Rq 
for dust that left the nucleus 0.5 day before perihelion to 
2.21 Rq for dust that left 0.3 day before perihelion and 
to 1.81 Rq for dust that left 0.1 day before perihelion. 
And, as expected, it is equal to the perihelion distance of 
the comet, 1.19 Rq, for dust that left at perihelion. Since 
the phase angle for particles on the syndyname /3 = 0.6 
released at perihelion is in Fig. 8 equal to 88°, nearly 
identical with the phase angle for particles released 0.7 
day before perihelion (see above), the utter absence of 
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Figure 10. Bizarre appearance of comet C/2011 W3 and its dust tail in an image taken with the COR2 coronagraph on board the 
STEREO-B spacecraft on December 16.517 UT, or 0.505 day after perihelion. The tail has two branches located, respectively, to the 
southwest and northeast of the Sun, while the comet's head is to the west-southwest of the Sun and north of the brightest part of the 
tail's southwestern branch. Careful inspection shows that the head's image is elongated along the line with the Sun, evidence that a new, 
postperihelion tail is being developed. The entire tail is fitted with a syndyname of /9 = 0.6, as shown in the panel that also provides the 
scale and orientation of the image. Only the heavily-drawn segments of the syndyname are seen as the tail in the image. The syndyname 
is calibrated by the times of dust release, reckoned in days from the time of the comet's perihelion passage (a negative number means 
before perihelion, and vice versa). Brightness-wise, the southwestern branch dominates because of forward scattering of sunlight. A 
broad secondary maximum is noted on the northeastern branch at a location corresponding to a release time of about 0.3 day before 
perihelion (December 15.7 UT). The tail then rapidly fades and disappears for a release time near 0.1 day before perihelion, suggesting an 
attentuation of activity hours before perihelion. Short segments of the syndynames f} = 0.4, 0.5, and 0.8 are also plotted, marked as A, 
B, and C, respectively, to illustrate a general deficit of dust with accelerations different from 0.6. The dots on these syndyname segments 
show the locations of particles released 0.1 day before perihelion. The northeastern tail's boundary is sharper to the south than to the 
north, indicating that dust with /3 2> 0.6 is more scarce than that with /3 < 0.6. No tail is apparent along the wide swath of space referring 
to release times during the first hours after perihelion. (Image credit: NASA/SECCHI consortium.) 



dust from emissions at perihelion has nothing to do with 
particle sunlight scattering effects. And since Figs. 8-10 
consistently show the absence of a tail for release times 
between approximately 0.1 day before perihelion and 0.1 
day after perihelion, that is, for particles whose perihe- 
lion distances were always smaller than ^1.8 Rq, the 
most likely explanation for this "activity attenuation" 
is the effect of sublimation of microscopic dust near the 
Sun, as discussed in detail in the following section. 

Probably the most diagnostic evidence for this effect 
is offered by the northeastern branch of the tail in Fig. 
10, which is made up of the dust that left the nucleus 
later than 2.2 days before perihelion. At the time of ob- 
servation, these particles were more than 3 Rq from the 
Sun; the phase angle reached 110° for those at the edge 
of the occulting disk but dropped rapidly along the tail 
to 70° about 1° from the Sun and was nearly constant 
around 50° along much of the rest of it. The apparent 
brightening at the location of particles released ^0.3 day 
before perihelion is thus likely to reflect their increased 
production around that time. Table 7 indicates that the 
dust at the edge of the occulting disk reached perihelion, 
at 2.7 Rq from the Sun, only about 0.45 day after the 



nucleus, while the dust at the far end of this tail's branch 
was nearest the Sun, at merely ~1.8 Rq, practically si- 
multaneously with the nucleus. This means that at the 
time of observation all particles along the northeastern 
tail were already moving away from the Sun. The fading 
and disappearance of this tail at a location populated by 
submicron-sized particles released ^0.08-0.1 day before 
perihelion and moving in highly hyperbolic orbits (eccen- 
tricity ^2.3 or higher) strongly indicate that their peri- 
helion distance of ~1.8 Rq is indeed the limit at which 
dust begins to sublimate profusely. 

The syndyname of /3 = 0.6 continues to extend in Fig. 
10 a little farther from the Sun than the disappearing 
northeastern tail and then turns sharply back, almost 
180°, running nearly parallel to the tail, in part behind 
the occulting disk, all the way to the comet's head. No 
trace of dust debris, all of which moved in orbits with 
perihelion distances smaller than ^1.8 Rq, is — just as 
in Figs. 8 and 9 — detected along this arc of the syn- 
dyname. Thus, no particles released between about 0.1 
day before perihelion and at least 0.1 day after perihelion 
appear to have survived, again pointing to ^1.8 Rq as 
the sublimation cutoff. 
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Table 7 

Orbital Elements of Dust Particles in the Tail of Comet C/2011 W3 



Orbital elements of released dust particles* as function of parameter f} 

Time Distance 

of from Perihelion distance (Rq ) Orbit eccentricity^ Perihelion time (days'*) 

release Sun 

(days'") {Rq) 0.3 0.6 1.0 1.5 2.5 0.3 0.6 1.5 2.5 0.3 0.6 1.0 1.5 2.5 



-10 


108.51 


1.69 


2.87 11.40 38.50 


-5 


67.94 


1.69 


2.81 


9.01 24.86 


-2 


36.37 


1.67 


2.69 


6.59 14.18 


-1 


22.51 


1.65 


2.55 


5.19 9.41 


-0.5 


13.80 


1.62 


2.37 


4.06 6.34 


-0.2 


7.09 


1.56 


2.07 


2.91 3.83 


-0.1 


4.22 


1.48 


1.81 


2.24 2.67 





1.19 


1.19 


1.19 


1.19 1.19 


+0.1'' 


4.22 


1.48 


1.81 


2.24 2.67 


+0.3 


9.55 


1.59 


2.21 


3.38 4.77 



65.98 1.0132 1.0790 1.1240 1.0241 

41.58 1.0212 1.1237 1.1921 1.0383 

22.60 1.0393 1.2212 1.3368 1.0704 

14.26 1.0628 1.3396 1.5074 1.1116 

9.01 1.1007 1.5154 1.7536 1.1766 

4.94 1.1884 1.8778 2.2456 1.3224 

3.17 1.3009 2.2894 2.7905 1.5023 

1.19 1.8570 3.9998 4.9999 2.3333 

3.17 1.3009 2.2894 2.7905 1.5023 

6.44 1.1428 1.6955 2.0001 1.2473 



+0.748 +1.862 +4.686 +4.148 -0.038 

+0.387 +0.967 +2.248 +1.931 -0.095 

+0.164 +0.407 +0.816 +0.653 -0.106 

+0.087 +0.206 +0.356 +0.259 -0.095 

+0.044 +0.098 +0.139 +0.082 -0.078 

+0.015 +0.027 +0.024 -0.002 -0.056 

+0.004 +0.003 -0.004 -0.018 -0.040 

0.000 0.000 0.000 0.000 0.000 

-0.004 -0.003 +0.004 +0.018 +0.040 

-0.025 -0.051 -0.060 -0.022 +0.066 



f' Release (ejection) velocity assumed to be zero. 

Reckoned from the time of perihelion passage of the comet; minus sign means before perihelion, and vice versa. 
*^ For /3 = 1 (motion in a straight line), the eccentricity is by definition infinitely large regardless of the release time. 

This entry is included to illustrate, by comparison with the entry —0.1 day, the symmetry with respect to perihelion. 



The panel in Fig. 10 also displays segments of the syn- 
dynames (3 = 0.4, 0.5, and 0.8 for a part of the north- 
eastern tail, where they do no overlap with the syndy- 
name j3 = 0.6. A clearly apparent property of this tail 
branch is its more gradual fading to the north of its sharp 
edge than to the south, where the drop is abrupt. Dust 
with (3 < 0.6, located to the north of the edge, has the 
perihelion distances systematically smaller, in agreement 
with the more limited presence of surviving dust than at 
/3 = 0.6. The tail no longer appears to contain particles 
with the acceleration parameter f3 significantly exceeding 
0.6, which should be located to the south. 

On the other hand, the tail's bright southwestern 
branch in Fig. 10 was made up of dust that left the nu- 
cleus between 8 and 3.6 days before perihelion. At the 
time of observation these particles were between 16 Rq 
and 5 Rq from the Sun. Table 7 shows that they passed 
through perihelion between 1.5 and 0.7 days after the 
comet, so they were all still approaching the Sun in Fig. 
10. And, as already alluded to above, this southwestern 
branch of the tail looks bright because of effects of for- 
ward scattering of sunlight, the phase angle varying from 
93° for dust released 8 days before perihelion to 134° for 
dust released 3.6 days before perihelion. At perihelion, 
the distance from the Sun was from 2.77 Rq to 2.86 Rq. 

The comet head's elongation in Fig. 10, already men- 
tioned above, indicates that dust particles were released 
from the nucleus starting about 4 hours, or slightly less 
than 0.2 day, after perihelion, and were therefore -^8 
hours old at the time of observation. This new, post- 
perihelion tail should have indeed extended from the nu- 
cleus in the antisolar direction. If consisting of particles 
with i3 = 0.6, its length at the time the image was taken 
should have been between 4' and 6'. 

There is no evidence in Figs. 8-10 for dust with /3 > 1. 
The question of what happened to it is again answered 
with the help of Table 7, which conveys two important 
facts on the whereabouts of particles in convex hyper- 
bolic orbits: (i) those released from the comet up to 
about 2 days before perihelion moved in orbits with per- 
ihelion distances always exceeding (often considerably) 



9 Rq , whereas (ii) more recent ones passed through per- 
ihelion nearly simultaneously with the comet, so that 
at the times the images in Figs. 8-10 were taken this 
dust already was at heliocentric distances much larger 
than the perihelion distance. In summary, dust particles 
with /3 > 1 did not contribute to the tail's postperihelion 
brightness because of their dispersal over a large volume 
of space and reduced light scattering efficiency far from 
the Sun. 

To address the issue of what kind of dust was released 
by comet C/2011 W3 after perihelion, we collected infor- 
mation on the tail length in the period of time from late 
December 2011 to mid-March 2012. A selection of 54 re- 
ported photographic and visual observations is listed in 
Table 8. Nearly all tail lengths from times after Decem- 
ber 21 can be explained by submicron-sized debris from 
the event of December 17.6 UT (Sec. 5), if the radiation- 
pressure accelerations /3 of up to 2.5 are allowed. On the 
other hand, constraining /3 to 0.6 requires that in most 
cases the tail (or at least its far reaches) derive from the 
activity prior to the event of December 17. To distin- 
guish between the two scenarios, it will be necessary to 
measure accurately the position angles of the tail, which 
after December 26 became almost perfectly straight. At 
present we are unaware of any such measurements. 

At the other end of the particle-size spectrum, a ques- 
tion arises as to how large the dust debris from the nu- 
cleus disintegration event of December 17 must be in 
order that it be released into elliptical orbits and not be 
lost to interstellar space. This can readily be determined 
by equating the orbital velocities of the comet and the 
debris on December 17.6 UT. For the debris to move be- 
low the escape limit, one requires that its /3cU be less 
than /3par for the parabolic limit, which follows from a 
condition 

— -(-) =(l-/3par)— , (18) 

''rel V^/ comet ''rel 

where (l/a)comct is the inverse semi-major axis of the 
comet's orbit and Vrei is the heliocentric distance at the 
time of release, rrei = 0.144 AU (Sec. 5). Thus, for the 
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Table 8 



Selection of Reported Postperihelion Tail Lengths of C/2011 W3 and Lengths Derived for Two Parameters /3o 



Observation 
time 
(Ui) 


Reported 
length 
of tail 


Predicted length^ 
/3o = 2.5 /3o = 0.6 


Required release time'' (UT) 
/3o = 2.5 /3o = 0.6 


Observer*^ 


Reference 


2011 Dec. 21.07 


10° 


o 

7.2 


o 

1.8 


Dec. 17.3 


Dec. 16.4 


N. Wakefield 


Mailing list 


21.30 


>5 


8.1 


2.0 




<Dec. 16.8 


*W. Souza 


Green (2012a) 


21.7 


13 


9.8 


2.4 


Dec. 17.3 


Dec. 16.4 


L. Barnes 


Mailing list 


21.73 


14.3 


10.0 


2.4 


Dec. 17.25 


Dec. 16.36 


R. McNaught 






lU 


12.7 


3.1 




Dec. 16.6 


W . oouza 


tureen (^zuizaj 


OO Qfl 


lo 


10 7 
IZ. / 


o. 1 


T^an 1 7 ^7 

uec 1 / .0 / 


uec iD.o 


*A. Aniorim 




22.32 


20 


12.8 


3.1 


Dec. 17.15 


Dec. 16.3 


*M.Goiato 


" 


22.7 


16 


14.8 


3.6 


Dec. 17.5 


Dec. 16.4 


L. Barnes 


Mailing list 


22.7 


15 


14.8 


3.6 


Dec. 17.58 


Dec. 16.5 


*D. Seargent 




23.31 


15 


18.3 


4.5 




Dec. 16.55 


*M. Goiato 


Green (2012a) 


ZO. 1 


OO 


20.7 


5.1 


Dec. 17.5 


Dec. 16.4 


J. iiiorooK 


Mailing list 




01 7 
zl. / 


on s 


0. 1 


uec 1 / .OO 


uec. ID. 4 


H. iviciNaugnt 




24.32 


20 


24.8 


6.1 




Dec. 16.5 


*M. Goiato 


Green (2012a) 


24.69 


28 


27.4 


6.8 


Dec. 17.57 


Dec. 16.4 


*D. Seargent 


Mailing list 


24.7 


27 


27.5 


6.8 




Dec. 16.4 


T. Barry 


" 


25.7 


30 


34.7 


8.8 




Dec. 16.5 


J. Dunphy 




26.60 


33 


41.2 


10.8 




Dec. 16.5 


* R. Kaufman 




26.63 


38 


41.^ 


1 n Q 




uec. ID. 4 


J. Drummond 




26.68 


37 


41.8 


11.0 




Dec. 16.45 


*D. Seargent 




26.71 


38 


42.0 


11.1 




Dec. 16.4 


M. Mattiazzo 




28.28 


30 


52.3 


15.1 




Dec. 16.8 


*W. Souza 




28.74 


30 


55.0 


16.4 




Dec. 16.85 


*A. Pearce 




on tin 


QO 


59.2 


18.8 




Dec. 16.9 


R. Kaufman 




on 


QO n 


Q 

oy.y 


1 Q 

±y . z 




uec iD.y 


*A. Pearce 




30.71 


25 


63.7 


22.1 




Dec. 17.4 


M. Mattiazzo 


" 


31.58 


30 


66.3 


24.7 




Dec. 17.3 


R. Kaufman 


" 


31.67 


45 


66.6 


24.9 




Dec. 16.75 


*D. Seargent 




31.72 


25 


66.7 


25.1 






*R. McNaught 




Om O Tr^Ti 1 TQ 

ZxjlZ Jan. i. /o 


OO 
ZZ 


68.8 


28.0 






rt. iviclNaugnt 




9 97 




Dy.D 


OQ A 

zy.4 






iVl. VjOld/TiO 


i^reen t^zvjizaj 


2.6 


>40 


69.9 


30.3 




<Dec. 17.1 


L. Barnes 


Mailing hst 


2.73 


36 


70.1 


30.6 




Dec. 17.3 


*R. McNaught 




3.73 


30 


70.7 


33.0 






*R. McNaught 




4.63 


33 


70.9 


34.8 






R. Kaufman 




A an 

4.oy 


oU 


70.9 


35.0 






IVl. iviattiazzo 




4. / U 


QQ 
OO 


1 u.y 


OO.U 




uec. 1 ( .4 


J. llIDrOOK 




5.70 


35.5 


70.6 


36.7 






J. Tilbrook 


" 


6.74 


21 


69.9 


38.0 






R. McNaught 


" 


13.48 


>18 


61.1 


39.2 






R. Kaufman 




14.48 


46 


59.7 


38.7 




Dec. 17.0 


R. Kaufman 




10. 4o 


>40 


58.2 


38.1 




<Dec. 17.0 


R. Kaufman 




10. 


A 1 

4 i 


58.2 


38.1 




Dec. 16.9 


L. Barnes 




16.48 


45 


56.8 


37.4 




Dec. 17.0 


L. Barnes 


" 


18.50 


37 


54.1 


36.0 




Dec. 17.5 


L. Barnes 


" 


20.48 


39 


51.6 


34.7 




Dec. 17.2 


R. Kaufman 


I 


22.5 


29 


49.4 


33.4 










23.5 


39 


48.4 


32.7 




Dec. 17.0 


L. Barnes 




24.5 


30 


47.5 


32.1 






L. Barnes 




25.55 


26 


46.6 


31.5 






R. Kaufman 




26.52 


31 


45.7 


31.0 






L. Barnes 




28.52 


17.5 


44.2 


30.0 






L. Barnes 




30.59 


26 


42.8 


29.0 






R. Kaufman 




Feb. 12.46 


7 


36.7 


24.8 






R. Kaufman 




Mar. 16 


1.2 


28.3 


18.9 






L. Barnes 


Barnes (2012) 



^ Assuming the tail consists only of dust released during the outburst that peaked on Dec. 17.6 UT, about 1.6 days after perihelion. 
^ For reference, the passage through perihelion occurred on Dec. 16.01 UT. 
^ Asterisk preceding the observer's name indicates a visual detection. 

^ Mailing list — ^http: / / tech.groups.yahoo.com/group / comets-ml|| messages 18968—19383. 
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debris in elliptical orbits 



1 

Pell < o^rcl - 

2 V a 



= 0.000915, 



(19) 



and the particle diameters Xdi at an assumed bulk den- 
sity of 0.4 g cm""^ are >0.3 cm, if released during the 
December 17 event. The correction to the barycenter of 
the solar system would change this result imperceptibly. 
The same condition for particles released at perihelion 
would yield /3cii < 0.000035 and x^w > 8.2 cm. We con- 
clude that only coarse-grain and larger debris continues 
to orbit the Sun in elliptical orbits. This is a way to gen- 
erate SOHO-like minicomets that move in orbits very 
similar to that of the comet from which they originate, 
except for the orbital period. 



11. SUBLIMATION OF MICROSCOPIC DUST NEAR THE 

SUN 

The existence of the so-called dust-free zone around 
the Sun, a result of profuse sublimation of interplane- 
tary dust, has been extensively studied. The problem is 
well reviewed by, for example, Mann ct al. (2004), based 
in part on earlier works by Krivov, Kimura, & Mann 
(1998), Kimura, Ishimoto, & Mukai (1997), and others. 
More recently this topic was also discussed by Kama, 
Min, & Dominik (2009) in connection with the problem 
of an inner boundary of protoplanetary discs. Mann et 
al. (2004) list the radius of the sublimation zone for dif- 
ferent materials, as computed by various authors. The 
numbers of particular interest are: 1.5-4 Rq for quartz 
(compact), 4 Rq for glassy carbon (both compact and 
fluffy), 5-6.5 Rq for Mg-rich pyroxene (both modes), 
10 Rq and 13.5 15.5 Rq for Mg-rich olivine (compact 
and fluify, respectively), 14 Rq for astronomical silicate 
(compact), and 11-24 Rq for a compact iron sphere. In 
most cases, however, the particles were assumed to orbit 
the Sun in a nearly-circular path until complete attri- 
tion, while for comctary particles moving along strongly 
hyperbolic orbits the critical perihelion distance is ex- 
pected to be even slightly less. 

The mineralogical and morphological properties of sili- 
cate grains in sungrazing comets were studied extensively 
by Kimura et al. (2002), who used a large set of light 
curves of the SOHO Kreutz minicomets as the data base. 
The authors concluded that the dust tails contained ag- 
gregates of submicron-sized crystalline grains, but not 
amorphous grains. The sublimation of fluffy olivine ag- 
gregates was proposed by them to explain the downturn 
on the light curves near 11-12 Rq. An increasing outflow 
of pyroxene fluffy aggregates was suggested as the source 
for a secondary maximum on the light curves near 7 -Rq . 
Their hypothesis does not appear to assign any signifi- 
cance to the range of heliocentric distances between 1.5 
and 2 Rq. 

To examine the sublimation properties of microscopic 
dust surviving in hyperbolic orbits with these small per- 
ihelion distances, consider a spherical dust particle of a 
bulk density p and radius cij-d released from the nucleus 
at time t^ci- Let the particle, which moves in a field of 
reduced effective gravity, sublimate at a temporally vari- 
able rate asubi(t)i such that at the time of observation, 



tobs) the radius is reduced to asn ^ Orei, so that 

rtohs 

Ctfin = Orel — J 0,suh\{t) dt, 



(20) 



where the sublimation rate is defined as a positive quan- 
tity. The rate dgubi depends strongly on the particle's 
equilibrium temperature T{t) and can be expressed in 
terms of the mass sublimation rate per unit area, 2suh\ 
(in g cm~^ s~^): 



<isubl(2^) 



--subl 



(T) 



P 



(21) 



where ZsuU is a function of the sublimation pressure, 

Psuhl{T), 



--subl 



(T) 



V27r3fjr 



Psubl(T), 



(22) 



yit is the molar weight (in g mol ^), Jft is the gas constant 
(in erg K~^ mol~^), and 



psubi(r) = Aexp 



L 



(23) 



with A and L, the latent heat of sublimation, being con- 
stants. 

The position at time t of a particle, moving in a con- 
cave hyperbolic orbit of the perihelion distance q and 
eccentricity e, is given by a hyperbolic eccentric anomaly 
/, which is related to t by 



esinh/ - / : 



A;gra^(l-/3)3(e-l)i 



{t-tl 



(24) 



where fcgrav is the Gaussian gravitational constant, /3 < 1 
is again the radiation-pressure parameter from Eq. (1), 
and is the time of the particle's passage through peri- 
helion. The anomaly / can for any t — t%he computed by 
successive iterations of Eq. (24). If t is reckoned from t^r, 
the time of the comet's passage through perihelion, then 
in Eq. (24) one substitutes f — t* = t — — (t* — f^). 
Since the heliocentric distance at time t is equal to 



r(i) = 



e- 1 



(esinh / - 1), 



(25) 



we can insert quantities from Eqs. (21) to (25) into Eq. 
(20) to obtain the integral in the form 



/■tobs f-i 

I disubi dt = ^l 



exp 



L 



rdf, (26) 



where /rei and fobs are the hyperbolic anomalies at times 
trei and fobs) and 



^ = 



A/i2 



[27r5ft(l-;3)(e- 1)]- 



(27) 



For a given tohs, the time of release t^ei and the 
radiation-pressure parameter /3 are determined by fitting 
the imaged tail's morphology with a set of syndynames 
that describe the distribution of dust particles in the tail. 
These computations automatically provide the particles' 
perihelion distance q and eccentricity e. Now the magni- 
tude of the integrated sublimation effect depends on four 
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quantities — the sublimation-pressure constant A (re- 
lated to the entropy of the particle material), the latent 
heat (or enthalpy) of sublimation L, the molar weight 
H, and the bulk density p. The effect is also a func- 
tion of the variations in the equilibrium temperature T 
along the particles' hyperbolic orbits. In practice, the 
computations are difficult because this temperature is a 
complicated function of the optical, thermophysical, and 
morphological properties of the particle material, all of 
which vary with time and heliocentric distance. 

We limit ourselves to a brief examination of the likeli- 
hood that the dust of C/2011 W3 that begins to subli- 
mate profusely at ^^1.8 Rq consists of olivine-dominated 
silicate particles released from the nucleus, respectively, 
0.08 and 0.12 day before perihelion and subjected to solar 
radiation pressure of ^ = 0.6. The first, more conserva- 
tive case implies that q = 1.7294 Rq, e = 2.4513, and 

— = —0.0008 day, whereas the other case leads to 
q = 1.8841 Rq, e = 2.1678, and t*^ - U = +0.0078 day. 
For the temperature we use only a fairly crude approx- 
imation based on a black-body temperature and an as- 
sumed thermal radiative regime, 



T(r) 



280 K 



r 2 



(28) 



The sublimation rate and other properties of olivine 
depend on the ratio of magnesium to iron in its formula 

Mg2Fe2t^Si04, where < z < 2. Following Barthelmy tion from the r^'^ dependence very close to the Sun, as 



plot in Fig. 11 the dependence of L on A/i 2 p i for the en- 
tire set of olivine members. In their study, Kimura et al. 
(2002) used Nagahara, Mysen, & Kushiro's (1994) sub- 
limation parameters for forsterite, but the molar weight 
for olivine with z = 1.1, and the density of 3.3 g cm~^, so 
that in Fig. 11 their point would be located just 0.04 in 
log(A/i2p^i) to the right of forsterite. For reference and 
orientation we also plot in the figure: (i) the data point 
for "silicates", as used by Kimura, Ishimoto, & Mukai 
(1997) — A = 1.07 X 10" Pa, L = 113.98 kcal mor^, 
^ = 67.0 g mol^^, and p = 2.37 g cm^"*: (ii) the num- 
bers that the most recent CRC Handbook for Chemistry 
and Physics (Haynes 2011) provides for quartz (Si02): 
A = 1.38 X 10" Pa, L = 134.55 kcal mol'S n = 60.08 g 
mol^^, and p = 2.65 g cm~^; and (iii) an extrapolated re- 
lationship between A and L, as used by Sekanina (2003) 
for some metallic elements in his study of the light curves 
of SOHO sungrazers, with an average value of 4 for the 
ratio Y^//9, which varied roughly between 1 and 7. 

To the extent that we can — despite the approxima- 
tions — express some confidence from the meaning of the 
plot in Fig. 11, we conclude that silicates dominated by 
Mg-rich olivine could represent a plausible candidate for 
the material that made up dust particles that began to 
sublimate profusely in the tail of comet C/2011 W3 near 
a heliocentric distance of 1.8 Rq. 

An effect that has not been accounted for in this pa- 
per is the deviation of the radiation-pressure accelera- 



(2010), we adopt z = 1.6 for the empirical formula, 
p = 153.31 g mol~^ for the molar weight, and p = 3.32 
g cm^'^, an average mineralogical density, for the bulk 
density of submicron-sized olivine-dominated grains. 

We now require that the integrated sublimation effect 
computed from Eq. (26) amounts to 0.2 pm in the parti- 
cle radius, the average size of silicate grains near ^ = 0.6 
(e.g., Kimura, Ishimoto, & Mukai 1997; Kimura et al. 
2002) and search for the relationship between the two re- 
maining unknowns — the sublimation-pressure constant 
A and the heat of sublimation L — that satisfy this total 
loss in the particle radius between trci and iobs- The two 
values of trei thus lead to two sets of solutions, which are 
insensitive to the observation time, as long as the ob- 
servation takes place after the particle's passage through 
perihelion. The solutions are presented in Fig. 11 in a 
range of 100 to 150 kcal mol~^ for the heat of sublima- 
tion and compared with several sources of information 
on sublimation of silicates and other potentially relevant 
materials. Because Eqs. (26) and (27) show that the in- 
tegrated sublimation effect depends on A only through 
a function Kp^ p~^ , it is this function that is plotted in 
Fig. 11 against the heat of sublimation L. 

Nagahara, Mysen, & Kushiro (1994) published the 
sublimation-pressure constants for the two endmem- 
bers of the olivine solid solution system. Including 
the molar weight and density, the data are as follows: 
A = 6.72 X 10" Pa, L = 129.74 kcal mol'S p = 140.69 
g mol~^, and p = 3.27 g cm~^ for forsterite (Mg- 
endmembcr, for which z = 2); and A = 2.48 x 10"'^^ Pa, 
L = 119.94 kcal mol-\ p = 203.78 g mor 1, and p = 4.39 
g cm~^ for fayalite (Fe-endmember, for which z = 0). Be- 
cause the sublimation properties vary smoothly between 
the two endmembers (e.g., Gail & Sedlmayr 1999), we 



briefly mentioned in Sec. 6. The reason for this devia- 
tion is twofold: (i) very close to the Sun, the particle is 
subjected to the radiation coming only from that part of 
the photosphere that is above the particle's horizon, and 
(ii) the photosphere is closer to the particle by up to 1 
Rq than the Sun's center. The first fact decreases the 
magnitude of the radiation pressure relative to that from 
the standard inverse-square power law, while the latter 
one increases it. Combined, the second influence pre- 
vails, and the integration over the relevant region shows 
that at a distance r from the Sun's center the actual 
acceleration, exceeds the standard value, (i, by 



^=2 



1- 1 1- 



Rq 



(29) 



It is clear that (3*^/3 for large values of the ratio 
t/Rq, while /?*//? is equal to 1.01 for r = 5 Rq, 1.15 
for r — 1.5 Rq, and 1.29 for r = 1.2 Rq. For the bor- 
derline heliocentric distances involved in the sublimation 
problem this effect does not exceed 10 percent and is not 
significant enough to invalidate the results. 

12. RAPID BUILDUP OF THERMAL STRESS IN NUCLEUS' 
INTERIOR AND PROCESS OF CATACLYSMIC 
FRAGMENTATION 

From the information presented, the process of dis- 
integration of comet C/2011 W3 appears to have been 
followed for at least four days immediately after peri- 
helion. A drop in activity around December 15.7 UT, 
or 0.3 day before perihelion, may have been related to 
sublimation of dust (Sec. 11) and was not necessarily an 
early signature of the process of disintegration. On the 
other hand, the sudden and permanent loss of the nuclear 
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Figure 11. The relationship between the parameter Afi^ p ^ — a function of the sublimation-pressure constant, molar weight, and bulk 
density — and the latent heat of sublimation L. The uninterrupted thick curve is the set of solutions that refer to the onset of profuse 
sublimation of dust in the tail of comet C/2011 W3 at 1.73 Rq from the Sun. The dashed thick curve is the set of solutions that refer 
to the onset of profuse sublimation at 1.88 Rq- The two squares are the locations in the plot of the two endmembers of the olivine solid 
solution system, based on the sublimation-pressure constants by Nagahara, Mysen, & Kushiro (1994): forsterite (Fo; Mg-endmcmber) and 
fayalite (Fa; Fe-endmember). The dots on the Fo-Fa connecting line refer, from left to right, to the members of the olivine system, whose 
formula is, respectively, Mgi.(iFeo.4Si04, Mgi.2Fc().gSi04, Mg().xFei,2Si04, and Mg().4Fei,6Si04. For reference and orientation, we also 
plot the data points for "silicates," as used by Kimura, Ishimoto, & Mukai (1997) (KIM), and for quartz, Si02 (Haynes 2011), and an 
extrapolated relationship between A and L, as used by Sekanina (2003) for some metallic elements in his study of the light curves of SOHO 
sungraaers, with an average value of 4 for the ratio ^fiij p- The constants for olivine used by Kimura et al. (2002) are not shown, as they 
nearly coincide with those for forsterite. 



condensation observed first between December 19.4 and 
20.3 UT (Fig. 1) was clear evidence of terminal collapse. 
There has been no trace of new activity detected from 
December 20 on, indicating that the comet's progressive 
fading after this date only reflected the rate of dispersal 
in space of the dust cjecta released during the December 
16-20 period of activity. 

It is not absolutely clear whether in this critical pe- 
riod of time the activity proceeded more or less continu- 
ously or was dominated by discrete outbursts other than 
the terminal event (s). Either way, this activity must ac- 
count for the formation of the quasi-parabolic envelope, 
encircling the spine tail, and also visible in its develop- 
ing phase in the images from December 17 through 19. 
Although the particle-ejection velocities in the spine tail 
were much lower than in the envelope (Sec. 8), the out- 
burst that triggered the formation of the spine tail was 
obviously the most devastating for the comet in that it 
was associated with the cataclysmic fragmentation of the 
entire residual mass of the nucleus. The low velocities of 
mostly large-sized debris may explain why to an earth- 
based observer the process appeared to have taken at 
least another two days to complete. Simple calculations 
support this scenario. Indeed, released on December 17.6 
UT and expanding at -^20 m s~^, large dust in the coma 



would occupy a volume of about 6000 km, or 10", in 
diameter on December 19.4 UT, when the comet still 
displayed a nuclear condensation (Fig. 1), and nearly 
10,000 km, or 17", in diameter on December 20.3 UT, 
when the condensation was already gone. The diameter 
of the bright condensation in the Malargue image from 
December 19 is in fact about 30", which suggests a ma- 
jor contribution of dust from the previous activity, while 
the breadth of the nascent spine tail in the Malargue im- 
age from December 20 is just about the expected 17". 
The bright streamer in the December 19 image — the 
early appearance of the spine tail — can be traced to 
a distance of at least 6' from the nucleus; its end point 
indicates the radiation-pressure acceleration parameter 
of /? ~ 0.12, that is, the presence of dust particles larger 
than 10 /im in diameter, when released during the De- 
cember 17.6 outburst. 

It is important to realize that the timing of this event 
offers unequivocal evidence of a gradual deterioration 
of the comet's health and thereby demonstrates its ap- 
preciable resistance to the hostile environment in close 
proximity of the Sun. This conclusion clearly shows 
that the nucleus could not possibly have been held to- 
gether merely by self-gravity. The comet's significantly 
delayed response does not favor models that describe the 
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comctary nucleus as a strengthless or very poorly ce- 
mented rubblc-pile structure, which — by virtue of being 
unconsolidated debris — reacts to tidal forces promptly 
(e.g., Asphaug & Benz 1994, 1996; Solem 1994). Ev- 
idence from comet C/2011 W3 thus leads to the same 
conclusion as evidence from the close-up imaging of the 
nuclei of 81P/Wild and 9P/Tcmpcl, which were indepen- 
dently found incompatible with the rubble pile model by, 
respectively, Brownlee et al. (2004) and Thomas et al. 
(2007). 

In a very recent study, Gundlach et al. (2012) explain 
the survival of sungrazing comets like C/2011 W3 within 
the Roche limit of the Sun by the counterpressure due 
to virtually isotropic outgassing from the micleus, which 
more than compensates for the tidal disruption force. 
For C/2011 W3 such a scenario is, however, unnecessary, 
because its nuclear diameter at perihelion was almost cer- 
tainly less than 1 km and the tidal stress therefore less 
than ~10 Pa (as seen from Fig. 2 in Gundlach et al.'s pa- 
per), a benign effect for any comet except an essentially 
strengthless one. A major problem with Gundlach et al.'s 
scenario is that it is inconsistent with the extensively ob- 
served duphcity of the sungrazing comet C/1965 SI (e.g., 
Pohn 1965, Thackeray 1965, lannini 1966, Lourens 1966), 
whose estimated diameter of less than 9 km (Knight et 
al. 2010) and perihelion distance of 1.67 Rq place the 
comet deep inside Gundlach et al.'s "safe" zone. 

What property was it then that made C/2011 W3 sur- 
vive perihelion on the one hand but suddenly disintegrate 
nearly 2 days later on the other hand? Always suspected 
in any siibstantially delayed response are the thermal ef- 
fects in the comet's nuclear interior, as transport of heat 
takes time. The resulting temperature changes, dramati- 
cally enhanced in the sungrazers near perihelion, impose 
a severe burden in the form of thermal stress, which may 
not only dwarf the tidal stress but could last for extended 
periods of time. 

To investigate heat transport effects in the nucleus of 
C/2011 W3, we have followed the procedure outlined and 
applied by Sekanina (2009) to study the giant explosions 
of 17P/Holmes. The approach assumes that the solar 
radiation impinging on the nucleus is spent — in rather 
extensive bare areas of the surface, devoid of ice — only 
on thermal reradiation and conduction of heat into the 
interior. This case is therefore, in principle, contrary to 
that considered by Gundlach et al. (2012). Being in- 
terested merely in basic, order-of-magnitude information 
on thermal stress, we have applied only the standard 
(isothermal) version of the one-dimensional heat-transfer 
problem in a spherical object. We have assumed that 
an initial central temperature of the nucleus is 60 K and 
a coefficient of effective thermal conductivity Kes = 0.2 
W m~^ K^^, which was found the most likely value to 
fit the recurrence of the giant explosions of comet 17P 
once every 115 years or so (Sekanina 2009) and which is 
also well within the interval compatible with the range of 
the thermal inertia that Davidsson, Gutierrez, & Rick- 
man (2009) derived in their analysis of near-infrared ther- 
mal emission spectra of features on the nucleus of comet 
9P/Tempel. With the numbers used for the specific heat 
and the bulk density, this value of K^s is equivalent to 
an effective thermal diffusivity of KeS — 0.006 cm^ s^^. 

The results of applying the heat-transfer equation to 
the nucleus of comet C/2011 W3 are presented in Figs. 



12 and 13. Figure 12, a plot of temperature against 
depth beneath the surface, exhibits an enormous peri- 
helion asymmetry. The plot shows for example that the 
same temperature reached 0.6 day before perihelion at 
a depth of less than 10 meters is reached 1.6 days af- 
ter perihelion at a depth of more than 40 meters. A 
significant conclusion from Fig. 13, a plot of the postper- 
ihelion temperature with heliocentric distance, is that 
exccipt in the topmost 15-meter-deep layer, the tempera- 
ture keeps increasing as the comet recedes from the Sun; 
it is still increasing at large heliocentric distances. In this 
regime, the effects of thermal stress must keep increasing 
and must further be enhanced by sudden activity from 
scattered reservoirs or "pockets" of highly volatile sub- 
stances, such as amorphous water ice, reached by the 
rapidly penetrating thermal wave deep inside the nu- 
cleus. Sooner or later, the combination of thermal stress 
and such explosive events is bound to have catastrophic 
consequences for the nucleus. We thereby find an en- 
tirely plausible mechanism for cascading fragmentation 
of sungrazers, large or small, well away from the Sun. 

The general equations for the radial and tangential 
components of thermal stress, cr||, a±, at distance h from 
the center of a sphere of radius h-k are given by (Timo- 
shenko & Goodier 1970) 
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is the thermal stress parameter and 



T(s) s^ds, 



(30) 



(31) 



(32) 



with T(s) being the temperature at distance s from the 
center. The object's thermal and mechanical properties 
determine the coefficient of linear thermal expansion oil i 
Young's modulus Ey, and Poisson's ratio v. From the 
first of Eqs. (30) we note that at the surface {h = h-k) the 
radial component cr|| of thermal stress is always nil. On 
the other hand, at the center the radial and tangential 
components have the same magnitude equal to 



lim fTii = lim a±^ 



2X 



hi 



(33) 



where Tq = T(0) is the central temperature. In addition, 
there exists a mean interior temperature (T) such that 



phi, i-hi, 

/ T{s)s^ds = {T) s'^ds= ^{T)hl, 
Jo Jo 



so that 



lim CT|| = lim a± = fx [{T) - To] . 
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As the distribution of temperature in the nucleus' inte- 
rior keeps steadily increasing with time and heliocentric 
distance following the perihelion passage, so do both (T) 
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Figure 12. Temperature variations with depth beneath the sur- 
face derived for the nucleus of comet C/2011 W3 with a maximum 
diameter of 600 meters. The one-dimensional heat-transfer equa- 
tion has been solved, averaging the incident solar energy over the 
surface, assuming that the energy is spent only on the conduc- 
tion into the nucleus and thermal rcradiation (bare surface, devoid 
of ices or other substances that would sublimate), and adopting a 
Bond albedo of 0.03, an cmissivity of unity, a coefficient of effective 
thermal conductivity of 2 X 10^ erg cm^^ s^-' K^-', a specific heat 
capacity of 8 x 10® erg g^^ K^^, and a bulk density of 0.4 g cm"''. 
The individual curves are identified by the time from perihelion: 
0.6 day before perihelion, at perihelion, and 0.1 and 1.6 days after 
perihelion, when comet C/2011 W3 was at heliocentric distances of, 
respectively, 0.0731, 0.00555, 0.0196, and 0.144 AU. Also marked 
are the crystallization temperature of amorphous water ice and the 
melting points of forsterite and fayalite. 
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Figure 13. Temperature variations with heliocentric distance at 
the surfa<;e and at 10 depths beneath the surface derived for the 
nucleus of comet C/2011 W3. The computations are based on the 
model and assumptions described in the caption to Fig. 12. The 
surface-temperature variations are shown for preperihelion (bro- 
ken) and postperihelion branches of the orbit, all other curves are 
postperihclion. The broken line parallel to the axis of ordinates 
shows the heliocentric distance 0.144 AU, at which the cataxjlysmic 
fragmentation event occurred on December 17.6 UT. 



(much more so than Tq) and the thermal stress in the 
center. The tendency toward a sungrazcr's breakup far 
from the Sun due to thermal forces is confirmed even 
without the existence of subsurface reservoirs of highly 
volatile ices. 

In Eqs. (30) we have divided the entire range of dis- 
tances from the center, < ft. < /i*, into n separate in- 
tervals hi-i < h < hi {i = 1, . . . ,n), where ho ^ and 
hn = hir, and fitted the temperature T from Fig. 12 in- 
side each interval by a polynomial of power m, 

m 

T{h) = Y,ak,ih^, i = l,...,n. (36) 

/c=0 

The integral $(/i^) is then given by 

n m 
i=l fe=0 



and for distance h from the center in a (j-l-l)-st interval, 
hj < h < hj+i (0 < j < n — 1), the integral ^{h) is equal 
to 

j TO 

m 

+i:jfii'^''^'-^^n- (38) 

fc=0 

Because Fig. 12 shows that the temperature changes very 
little at depths exceeding 60-70 meters, the described 
approach can also be used to approximate the magnitude 
and distribution of thermal stress in the interior of the 
nucleus of smaller dimensions, whose center is assumed 
to be at depth /i^ — h rather than at hi^. Reckoning 
distance h (but not the interval boundaries hi-i and hi) 
from the new center at this depth, we can write Eq. (36) 
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in the form 



^C*) = + ^)'' =Y1 ^'''i^^' i=j,...,n, 



k=0 



fc=0 
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when h is in an interval hj-i < H < hj (1 < j < n), and 
for < fc < m 
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CoefBcients bk,i now replace ak,i in Eqs. (37) and (38), 
in which the interval boundaries have to be modified ac- 
cordingly. 

This set of formulas can also be used when the size 
of the nucleus is truncated to h» < /i^ on the outside, 
that is, when it is desirable to eliminate the contribution 

to thermal stress from a surface layer once it has been 
removed by erosion or otherwise destroyed. In this case 
the integral $(/i») is computed from Eq. (38) rather than 
(37), either with the coefficients ak,i or bk,i, depending 
on whether H is or is not nil. 

To allow for a gradual erosion of the surface layer in 
the proximity of the Sun, the internal-temperature dis- 
tribution presented in Figs. 12 and 13 has been used to 
compute the thermal-stress profiles for the nucleus 400 
meters in diameter at 0.6 day before perihelion, 280 me- 
ters in diameter at perihelion, and 150 meters in diam- 
eter 1.6 days after perihelion. The results are presented 
in two steps. In the first step, the relative variations in 
the radial and tangential components of thermal stress, 
normalized to its magnitude at the center of the nucleus 
at perihelion, are displayed against depth beneath the 
surface in Fig. 14. The major properties of the stress 
curves are as follows: 

(1) Just as with temperature, thermal stress continues 
to climb even as the comet recedes from the Sun; mea- 
sured by its magnitude at the nucleus' center, the ther- 
mal stress at perihelion is about 16 times higher than 0.6 
day before perihelion, but 2.4 times lower than 1.6 days 
after perihelion. 

(2) At a given time, the radial component of thermal 
stress is always tension, increasing from zero at the sur- 
face, as already noted; the tangential component reaches 
its maximum at the surface as compression, but changes 
to tension a few tens of meters beneath the surface, in- 
creasing toward the center. 

(3) Both components reach maximum tension at the 
center of the nucleus, where they have the same magni- 
tude, as implied by Eq. (33). 

(4) Cracks, whose formation is triggered by increasing 
thermal stress at large depths due to the thermal wave's 
gradual penetration into the interior, are bound to open 
preexisting subsurface reservoirs of highly volatile ices 
(such as amorphous water ice), initiating their explosion 
when heated up to the activation temperature. 

(5) The steady increase in thermal stress at the cen- 
ter of the nucleus over long periods of time after perihe- 
lion is likely to be a driving force behind the episodes of 
seemingly spontaneous cascading fragmentation at large 
heliocentric distances, which account for the SOHO sun- 
grazers' observed scatter in elements other than orbital 
period. 
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Figure 14. Relative thermal stress in the interior of comet C/2011 
W3 as a function of depth beneath the surfa<;e. Derived from the 
temperature variations 0.6 day before perihelion for the nucleus 
400 meters in diameter (solid), at perihelion for 280 meters in di- 
ameter (broken), and 1.6 days after perihelion for 150 meters in 
diameter (solid), the stress curves are normalized to thermal stress 
at the nucleus' center at perihelion. The segments of the curves 
are marked to distinguish between the radial and tangential stress 
components and whether the stress was a compression or tension. 



Measured in absolute units, the magnitude of thermal 
stress depends on the stress parameter x, given by Eq. 
(31) and consisting of three physical quantities: Young's 
modulus Ey, the coefficient of linear thermal expansion 
aL, and Poisson's ratio v. Since their values for cometary 
material are unknown, wc make in the following an effort 
to constrain each of them to our best ability by using 
metcoritic, lunar, and terrestrial analogues. Given our 
conclusion in Sec. 11 on the nature of dust sublimating 
from comet C/2011 W3 near the Sun, a plausible can- 
didate among the terrestrial samples appears to be an 
olivine-based silicate material. 

Young's modulus. There is a large number of papers 
with information on Young's modulus, but they mostly 
refer to samples of zero or near-zero porosity. For the 
Moon, Pritchard & Stevenson (2000) estimated that the 
effective Young's modulus in the outer layers of the lu- 
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nar lithosphcrc is between 5 and 10 GPa. For ordinary 
chondrites the topic was recently updated by Kimber- 
ley & Ramesh (2011). They also provided the results 
of their experiments aimed at determining the elastic 
and mechanical properties of a particular L5 ordinary 
chondrite found in Antarctica, showed that the outcome 
depends on whether the data were obtained under q\m- 
sistatic or dynamic conditions, compared their results 
to the findings from the numerous studies of the atmo- 
spheric breakup of meteors, and emphasized that the 
propagation of cracks is gradual because the crack speeds 
are much lower than the limit defined by elastodynamic 
fracture theory. Although Kimberley & Ramesh did not 
consider thermal stress, their conclusions are relevant, 
because they found that the value of the Young's mod- 
ulus for a particular sample with a 13 percent poros- 
ity depended on a strain rate and range of strain; at 
slow rates and strains of up to 0.01, Ey = 3.2 GPa, at 
strains ranging from 0.01 to 0.03 and much higher strain 
rates (appropriate in hypervelocity impact phenomena), 
Ey = 8.5 GPa; this behavior has also been observed in 
some brittle terrestrial materials, e.g., ceramics (Wang 
& Ramesh 2004). Pertinent to our problem is the lower 
value of Ey ■ 

From analysis of limited samples of meteorites, rarely 
including carbonaceous chondrites (e.g., Flynn, Klock, & 
Krompholz 1999), most values found for Young's modu- 
lus were one order of magnitude higher than the quoted 
result by Kimberley & Ramesh, but their porosities were 
mostly lower than 13 percent. A key work was pub- 
lished by Yomogida & Matsui (1983) on the mechani- 
cal, thermal, and elastic properties (measured quasistat- 
ically) of 20 ordinary chondrites of H and L types. For 
some of these properties the authors compared the me- 
teorites with lunar samples studied by others. One of 
Yomogida & Matsui's findings was the strong depen- 
dence of the elastic properties on porosity. They showed 
that even a low porosity reduces the values of the elastic 
parameters, including Young's modulus, much more for 
the meteorites and the lunar samples than for terrestrial 
rocks. While there are fairly large differences between the 
individual samples, the steep drop in Young's modulus 
Ey{'>Ij) with porosity ij} is for 19 out of the 20 chondrites 
fitted by 



Ey{iI>) = const -exp 



-U 



(for V > Vo), (41) 



where U = 13.8 ± 1.2 and const = 85^11 x lO^^ Pa 
{tpo = 0.05) for 13 samples (high Ey group) and 
X 10^2 Pa (V'o = 0.02) for 6 samples (low Ey 
group). This expression provides a siibstantially better 
fit than Warren's (1969) formula used by Yomogida & 
Matsui (1983) in their Fig. 8 to match a normalized 
elastic modulus of both the meteorites and the lunar 
samples. Although by 2001 the porosity (with values 
of up to 30 percent) had already been known for more 
than 450 stony meteorites (Britt & Consolmagno 2003), 
information on their elastic moduli seems to be lagging 
behind. 

A variety of studies of terrestrial samples provides con- 
straints on Young's modulus of olivine. However, many 
measurements refer to temperatures in excess of 1000 K. 



For example, Tait (1992) used a value of 197 GPa at 
1470 K, while Hiraga, Anderson, & Kohlstedt (1993) gave 
values of 159, 155, and 154 GPa at, respectively, 1373 K, 
1473 K, and 1523 K. Chung (1970) measured the bulk 
modulus and Poisson's ratio of olivine as a function of the 
Mg/Fe ratio at room temperature (296 K) and zero pres- 
sure. In terms of Young's modulus, his results show a sys- 
tematic trend from a maximum of 198 GPa for forsterite 
to a minimum of 141 GPa for fayalite. An extensive 
study of the dependence on the temperature of the elas- 
tic properties of a single crystal of forsterite was pub- 
lished by Suzuki, Anderson, & Sumino (1983), who found 
that the elastic moduli increase with decreasing temper- 
ature at a rate of about 1.3 percent per 100 K between 
OK and 1200 K. From their values of the bulk modulus 
and Poisson's ratio, Young's modulus of forsterite comes 
out to be 175 GPa at 1500 K and 207 GPa at 100 K, 
an increase of 18 percent. All these numbers are fairly 
consistent with the highest values of Young's modulus 
near 140 GPa measured by Yomogida & Matsui (1983) 
for the meteorites of low porosity. It therefore appears 
that the temperature effect is insignificant relative to the 
porosity effect. A detailed investigation of Young's mod- 
ulus variations with porosity for a variety of terrestrial 
materials was published by Phani & Sanyal (2005), who 
determined that for colloidal gel derived silica the rela- 
tive modulus varied as 



EyW _ (1 - 
(-E^y)o 



1 + iV 



(42) 



for porosities smaller than 0.85. The slope of the curve is 
somewhat similar to that shown by Yomogida & Matsui 
(1983) for a sintered material and tuff in their Fig. 8; it is 
considerably less steep than the slope for the meteorites 
and the lunar samples. 

Coefficient of lineEir tiiermal expansion. This coefla- 
cient is on a recently compiled "wish" list of 12 physical 
quantities of meteorites, the data on which are desired 
most by the scientific community (Consolmagno, Britt, 
& Opcil 2010), a fact that cogently illustrates the urgent 
need for this practically nonexistent information. 

It has been known for nearly a century that the thermal 
expansion coefficient is generally temperature dependent 
and that its rate of variation in crystallic solids is, at 
least at relatively low temperatures, approximately pro- 
portional to that of the specific heat capacity at constant 
volume (e.g., Austin 1952). Given that our primary in- 
terest is the temperature range from ~70K to ~300K, 
we have focused on a search for relevant sets of data on 
olivine-based silicates that should be as consistent with 
this temperature interval as possible. Once available, the 
data on the coefficient of linear thermal expansion have 
been fitted to satisfy a properly scaled formula for the 
Einstein crystal model (e.g., Rogers 2005), 



aL{T) = A 



sinh 



(43) 



where A and are constants to be determined from the 
data. (The known fact that this equation does not predict 
correct values of heat capacity at low temperatures is 
irrelevant to our using it to fit a set of cxl values.) 

Samples of a hot-pressed olivine aggregate were used 
by Aizawa, Ito, & Tatsumi (2001) to measure its vol- 
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umetric coefficient of thermal expansion, ay, but the 
temperature range was from 300 K up, which makes the 
data, requiring much extrapolation, only marginally use- 
ful. After conversion to aj, = ^ay, we found that fitting 
the data between BOOK and 600 K by Eq. (43) requires 
A = 12.0 X lO-'^K-i and O ^ 304 K. With these con- 
stants aL equals 1.02, 5.85, and 8.62 x lO-^K"! at the 
temperatures of, respectively, 100 K, 200 K, and 300 K. 

Another set of volumetric thermal expansion coeffi- 
cients is found in the paper by Suzuki, Anderson, & 
Sumino (1983), already mentioned in the subsection on 
Young's modulus. Their data refer to the crystal of 
forsterite and cover the desired range of temperatm-es. 
A fit by Eq. (43) requires that A = 9.14 x 10"^' K^^ and 
fl = 251.4K. The fit gives aL equal to 1.53, 5.54, and 
7.27 X 10"^ K-i at 100 K, 200 K, and 300 K, respectively. 
The fit to the measured values is perfect for the first two 
temperatures, but too low by 12 percent for the third. 

Unfortunately, most information available on the ther- 
mal expansion coefficient refers to geology of the earth's 
crust, volcanology, and related scientific fields (e.g., 
Afonso, Ranalli, & Fernandez 2005). As a result, one sel- 
dom finds data pertaining to low temperatures and low 
pressures. In addition, the effect of porosity on thermal 
expansion is seldom addressed, and when it is, the re- 
sults and implications depend on the type of material and 
physical conditions (e.g., Faivre, Bellet, & Dolino 2000; 
Moretti et al. 2005; Hunter & Brownell 2006; Ghabe- 
zloo & Sulem 2009; Ghabezloo 2010), although generally 
there appears to be a tendency for the thermal expansion 
coefficient to increase with increasing porosity. 

Poisson's ratio. Besides the elastic moduli, Yomogida 
& Matsui (1983) also measured Poisson's ratio v for the 
already mentioned 20 meteorites. With the exception of 
a single anomalous sample (the same one that did not fit 
the; depcindence of Yoimg's modiilus on porosity), all 19 
remaining had their values of i/ between 0.12 and 0.29. In 
this data set there is a slight tendency for u to decrease 
with increasing porosity, but the scatter is too large for 
a functional fit; on the average, (j^) = 0.21 ± 0.05. 

For olivine, Chung (1970) found that Poisson's ratio 
was increasing systematically with the Fe/Mg ratio from 
0.242 for forsterite to 0.308 for fayalite. Just as with the 
thermal expansion coefficient, the dependence of Pois- 
son's ratio on porosity for terrestrial materials appears 
to vary from case to case. Gel derived silica, for example, 
has Poisson's ratio nearly constant between 0.15 and 0.19 
at porosities smaller than 0.6, but it then rapidly climbs 
to 0.26 at porosity 0.82 (Ashkin, Haber, & Wachtman 
1990; Phani & Sanyal 2005). On the other hand, Pois- 
son's ratio of sintered iron decreases from 0.30 at zero 
porosity to 0.20 at porosity 0.2 (Kovacik 2005). 

The temperature dependence of Poisson's ratio for 
forsterite is shown by Suzuki, Anderson, & Sumino's 
(1983) data. The rate amounts to <0.001 per 100 K and 
can be neglected. Compared to Young's modulus and 
the thermal expansion coefficient, Poisson's ratio has a 
relatively minor effect in computing thermal stress. 

Adopting stress parameter x for C/2011 W3. 
Whereas the compiited relative variations in thermal 
stress in the interior of comet C/2011 W3 in Fig. 14 
suggest an enormous growth near and after perihelion, it 
remains unclear whether the effect is sufficient to crumble 



and collapse the nucleus. Investigations of this kind can 
only be conducted in a framework of a particular comet 
model, with realistic estimates for the quantities that en- 
ter the stress parameter x- The extensive discussion of 
Young's modulus, the thermal expansion coefficient, and 
Poisson's ratio was intended to facilitate this effort. 

To anchor our efforts in the framework of a numerical- 
simulation construct for comet formation, evolution, and 
morphological makeup, we focus on Lasue et al.'s (2009, 
2011) work on aggregation and coUisional interaction of 
cometesimals in the primordial solar nebula and the in- 
ternal structure of cometary nuclei. Their developed 
model is based on a homogeneity exponent determining 
an aggregation regime, accounts for disruptive, sticking, 
compaction, and sintering processes, allows comparison 
with observational constraints and future in situ observa- 
tions, and provides predictions for the cohesive strength 
and the radial porosity profile of the nucleus' interior. 
In our choice of parameters we conciir with Lasue et al., 
who clearly prefer the layered striicture of the nucleus, 
which is determined by the range of 0.4 to 0.6 for their 
homogeneity exponent and which offers the characteris- 
tic values of the porosity and cohesive strength that we 
use in this study. 

Keeping in mind the complete disintegration of C/2011 
W3 after perihelion, our primary interest is the temporal 
variation in thermal stress deep inside the nucleus, which 
we refer to as the central thermal stress. By this term we 
understand stress in the volume of the nucleus that does 
not include the surface layer, where the thermal stress 
profile is very different from that at greater depths. In 
Fig. 14 the boundary between the surface layer and the 
deep interior approximately coincides with the apparent 
minimum on the stress curve, associated with the tran- 
sition of the tangential component a± from compression 
to tension. 

The extensive areas of the surface and the adjacent, 
relatively shallow subsurface layer of the nucleus that 
are presumed to be devoid of ices, must in the proxim- 
ity of perihelion suffer from even higher levels of ther- 
mal stress than the rest of the nucleus (Fig. 14). They 
are primary candidates for thermal-stress damage, being 
probably riddled with crisscross cracks and the debris 
removed soon after perihelion. These areas may in fact 
provide for the comet's earliest postperihelion activity. 
We will not deal with them in this paper. 

In Fig. 14, a rather interesting property of the thermal 
stress distribution in deeper layers, say, from a depth of 
~30 meters on, is that the magnitude of the total effect, 

a = y^^f + ai, (44) 

is almost independent of depth, especially at perihelion 
and afterwards. Because a fairly rapidly converges to 
(T|[\/2 = o-±V^ at the nucleus' center, it is both conve- 
nient and appropriate to deal with this limit as a mea- 
sure of thermal stress in our consideration of its temporal 
variations in absolute units — hence the reason for the 
term central. 

The discussion of the quantities that make up the stress 
parameter x suggests that in an effort to find their most 
appropriate values one will encoimter considerable diffi- 
culties and uncertainties. The two most severe problems 
are the strong dependence of Young's modulus on poros- 



28 



Sekanina & Chodas 



ity and the thermal expansion coefficient's variation with 
temperature. In the former case, there is the need to 
extrapolate the data way outside the covered range of 
porosity. In the latter instance, it is troubling to find 
a parameter, which the thermal stress theory handles 
as a constant, to be so strongly temperature dependent. 
Only the selection of a value for Poisson's ratio appears 
to be somewhat less controversial, thanks to this quan- 
tity's relative invariability. Yet, for all three quantities it 
is necessary to resort to other classes of objects and/or 
substances to approximate cometary material. Also, we 
have no choice but to rely on a comet model to estimate 
the porosity in the interior. With these caveats in mind 
we proceed to the next step. 

We note that Lasue et al.'s (2011) aggregate model pre- 
dicts a porosity close to 0.65 at large depths and near the 
center of a cometary aggregate regardless of whether the 
size distribution function of the cometesimals is Gaus- 
sian or a power law. Equation (41) then yields Ey = 8.6 
MPa from the high Ey group of meteorites and Ey = 3.1 
MPa from the low Ey group. For the coefficient of linear 
thermal expansion we have to adopt a mean value 
from a range of less than 1 to at least 6 or 7 x 10~^ K~^, 
based on Suzuki, Anderson, & Sumino's (1983) data for 
a single crystal of forsterite. Because of the possibility 
that an explosion from a reservoir of amorphous water 
ice at depths in excess of 40 meters was a trigger of, or 
a contributor to triggering, the terminal outburst (Fig. 
12), we have chosen q:l = 3 x 10~^K~^, which is near 
the mean and is also attained at 130 K, the temperature 
of activation of the ice's exothermic process of crystalliza- 
tion (e.g., Schmitt et al. 1989). Because of the dearth of 
relevant information, we have not included any porosity 
correction. For Poisson's ratio we have adopted f = 0.1, 
which is based primarily on the sample of meteorites and 
allows for an assumed modest decrease with porosity. Its 
value has the role of a minor correction to the stress pa- 
rameter, which comes out to be x = 29 Pa and 10 Pa 
with the two values of Ey. 

To address the issue of whether the magnitude of ther- 
mal stress determined in this way for C/2011 W3 could 
lead to the crumbling and collapse of this comet's nu- 
cleus, we have compared our results with the expected 
cohesive strength of a cometary aggregate provided by 
Lasue et al.'s model. Specifically, we have taken approx- 
imate upper and lower limits near the homogeneity ex- 
ponent of 0.5, as presented in Fig. 7b of Lasue et al. 
(2009), because they alluded to this point only in gen- 
eral terms in their second paper (Lasue et al. 2011); the 
understanding the reader gets is that the introduction 
of a size distribution function of cometesimals had no 
significant eff'ect on the cohesive strength. 

The resulting thermal stress curves, referring to the 
two adopted values of Young's modulus, are plotted 
against time from perihelion in Fig. 15. In broad terms, 
we conclude that thermal stress developing in the interior 
of the nucleus is: (1) insufficient to disrupt the object be- 
fore perihelion; (2) comparable to the cohesive strength 
in the immediate proximity (within hours) of perihelion, 
thus weakening the nucleus' structure; and (3) greater 
than the strength of the comet after perihelion to the 
extent that it can no longer hold together. 

Because of the uncertainties involved, one needs to ex- 



ercise a caution and emphasize that this conclusion is 
model and parameter dependent. Nevertheless, in the 
case of the thermal expansion coefficient, its chosen con- 
stant value has a tendency to suppress the degree of 
asymmetry of thermal stress relative to perihelion, that 
is, it overestimates thermal stress before perihelion and 
underestimates it after perihelion. As a matter of fact, 
this relation can be quantified; from Eq. (43) it follows 
that a change AT in the temperature leads to the fol- 
lowing change Aql in the coefficient of linear thermal 
expansion: 









2 


p coth. 











This formula indicates that, for example, a modest tem- 
perature increase of 5 K entails an increase in the rela- 
tive rate of the aL coefficient and therefore in thermal 
stress — of 2.4 percent at 200 K, 5.3 percent at 150 K, 
15.5 percent at 100 K, and an enormous 37 percent at 
70 K. This effect comes on top of the thermal stress in- 
crease that we have discussed below Eq. (35), and en- 
hances the role of thermal stress as the driving force be- 
hind the events of cascading fragmentation at large he- 
liocentric distances, as mentioned in point (5) between 
Eqs. (40) and (41); it also aids the role of thermal stress 
in unlocking potential reservoirs of highly volatile ices in 
the nucleus' interior. 

In summary, our suspicion that the delayed response of 
the nucleus of C/2011 W3 to the extremely high temper- 
atures in the proximity of perihelion was due to heating 
its interior and to the resulting effect of thermal stress 
finds support in this comprehensive analysis, which is 
limited only by incomplete information on cometary ma- 
terial the problem that no hypothesis can avoid. If 
our scenario is valid at least with an order-of-magnitude 
accuracy, it strengthens the notion that C/2011 W3 was 
able to withstand thermal stresses on the order of sev- 
eral kPa and therefore possessed far more cohesion than 
necessary to avoid collapse and disintegration right at 
perihelion, in the inner solar corona. 

13. STANDING OF C/2011 W3 IN THE KREUTZ SYSTEM 
AND ITS POSSIBLE PAST EVOLUTIONARY PATH 

The place of the now defunct comet C/2011 W3 in 
the hierarchy of the Kreutz system of sungrazers cannot 
be pinpointed with certainty. However, within a given 
model of the system, the comet's evolutionary path can 
be traced in some detail. The exercise involves modeling 
and interpretation of long-term changes in (a) the or- 
bital period and (b) the angular elements and perihelion 
distance. As amply explained in Sec. 1, the first issue 
is addressed in terms of the tidal-assisted fragmentation 
of the early parent objects of C/2011 W3 in the close 
proximity of perihelion, while the second issue involves 
a sequence of fragmentation episodes of its more recent 
precursors far from the Sun. 

As an example, we elaborate on the evolutionary path 
of C/2011 W3 in the context of the hierarchical scenario 
B of the Kreutz system, introduced in our recent pa- 
per (Sekanina & Chodas 2007). This scenario is based 
on the postulated identity of the spectacular sungrazer 
X/1106 CI with a comet' of A.D. 467. The true orbital 
period of this progenitor (or superfragment, in the termi- 
nology used in Sekanina & Chodas 2004) was obviously 



Comet C/2011 W3: Physical and Orbital Study 



29 



- 1 I I I I I I I I I I I I I I I I I I I I I I I I I I I - 

: THERMAL STRESS NEAR CENTER : 
- OF NUCLEUS OF C/2011 W3 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I 
-0.5 +0.5 +1 +1.5 

TIME FROM PERIHELION (days) 



Figure 15. Computed thermal stress near the center of the nucleus of C/2011 W3. The solid curves show the stress growth through 
perihelion for the two adopted values of Young's modulus, the horizontal dashed lines mark the estimated limits on the cohesive strength of 
the nucleus (from Lasuc ct al. 2009, 2011). The stress curves represent a conservative estimate; because the thermal expansion coefficient 
increases with the temperature, the actual postperihelion growth of thermal stress in the comet's interior is even steeper. The time of the 
terminal outburst of C/2011 W3, 1.6 ± 0.2 days after perihelion, is shown by the solid vertical line, its uncertainty by the parallel broken 
lines. 



639 years. The integration of the orbit of C/2011 W3 
back in time indicates (Table 4) that the previous return 
to perihelion occurred in January 1329 (vi^ith a formal 
uncertainty of ±2 years, which for this exercise is not 
critical). It should be emphasized that the 1329 orbital 
set is valid only in the unlikely absence of any fragmen- 
tation events involving the comet's precursors during the 
entire revolution about the Sun. Only in this hypothet- 
ical case does the true orbital period come out to be 
almost exactly 683 years. In order for a fragment an 
early parent of C/2011 W3 — of comet 467 to arrive at 
perihelion in 1329, its true orbital period must have been 
862 years, or 223 years longer than that of the principal 
fragment, comet X/1106 CI. This increase in the orbital 
period requires that the parent fragment of C/2011 W3 
split off from the main body of comet 467 near perihe- 
lion (say, at a heliocentric distance of 0.01 AU) with a 
separation velocity of ~2.2 m s~^ essentially in the direc- 
tion of the orbital motion (Table 8 of Sekanina 2002a). 
When this parent fragment arrived at perihelion in 1329, 
its orbital period was (with some uncertainty due to the 
indirect planetary perturbations) still 862 years. Since 
the true orbital period of C/2011 W3 is 683 years, the 
necessary decrease of 179 years in the orbital period re- 
quires that a new fragment broke off from the parent 
fragment in 1329, again in the immediate proximity of 
perihelion, with a separation velocity of ^3.6 m s^^ es- 
sentially in the direction opposite the direction of orbital 
motion. Separation rates of fragments corresponding to 



relative velocities of up to 5 or 6 m s~^ are known to 
be necessary for explaining the tidal-assisted splitting of 
sungrazing comets (Sekanina & Chodas 2007), so the re- 
quired velocities of less than 4 m s^^ provide a plausible 
evolutionary model for comet C/2011 W3. Similar num- 
bers would result if the adopted scenario B should be 
replaced with another one. 

Because of potential differences in the angular elements 
and perihelion distance between the orbits of C/2011 W3 
and the new fragment referred to above, it is highly un- 
likely that the two objects are identical. Much more 
probable is that the new fragment was a precursor to 
C/2011 W3, which between 1329 and 2011 underwent 
additional, nontidal fragmentation at large heliocentric 
distances, and that during these episodes one of the frag- 
ments acquired the necessary orbital orientation and per- 
ihelion distance to become C/2011 W3. If the first pre- 
cursor's angular elements in the 14th century were in the 
range of most SOHO sungrazing comets, that is, close to 
those of C/1963 Rl, the differences relative to C/2011 
W3 may have been substantial, about 35° to 40° in the 
argument of perihelion and the longitude of the ascend- 
ing node, some 10° in the inclination, and up to 10 per- 
cent of the Sun's radius in the perihelion distance. To 
bridge the gaps of this magnitude requires several non- 
tidal fragmentation events in the general proximity of 
aphelion (at heliocentric distances of, typically, 100 AU 
or so) at an assumed average separation velocity of 2-3 m 
s~^, a requirement by no means excessive. A simulation 
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Figure 16. Possible evolutionary path for comet C/2011 W3 in a broader context of the model scenario B (Sekanina & Chodas 2007), 
which assumes that comet C/1843 Dl, rather than C/1882 Rl, is the most massive fragment of X/1106 CI. The parallel fragmentation 
branch involving C/1882 Rl is not shown in this schematic representation. The diamond-shaped boxes describe events of the tidal-assisted 
and cascading nontidal fragmentation processes, the rectangular boxes are their products. The separation velocities involved arc given as 
Ay; the sign refers to the direction along the orbital-velocity vector: positive means ahead, negative behind. The fragmentation hierarchy 
is presented schematically, with only some paths among the great multitude of possible chain events being depicted. 



model for a sequence of nontidal fragmentation events 
scaled for the Kreutz system from the case of D/1993 
F2 (Shoemaker-Levy 9) shows that the number of these 
episodes can easily exceed a dozen during one revolution 
about the Sun (Sekanina 2002a, 2002b). Their high fre- 
quency is also obvious from the arrival rate of the SOHO 
minicomets, the high-generation products of the cascad- 
ing fragmentation process, the age of each of them being 
manifestly shorter than one revolution about the Sun 
(e.g., Sekanina & Chodas 2004). 

In the broader context of the fragmentation hierarchy 
of the branch of the Kreutz system, which in our Scenario 
B also contains the sungrazcr C/1843 Dl, the described 
evolutionary path of C/2011 W3 is depicted in Fig. 16. 
In this more general scheme, the parent fragment may 
have reached perihelion years or even decades before or 
after 1329, because fragmentation at large heliocentric 
distances entails modest changes in the orbital period of 
up to a few years per event (Table 8 of Sekanina 2002a). 

Regardless of the details of the evolution of C/2011 
W3, its origin appears to be linked to the expected 
new, 21st century cluster of bright sungrazers. What 
remains unclear is its subgroup membership. While al- 
ready Kreutz (1901) was thinking in terms of two sub- 
groups when investigating historical sungrazers, this divi- 
sion for the bright members of the Kreutz system became 
a hit after Marsden's (1967) publication of his classical 



paper, which listed C/1843 Dl, C/1880 CI, and C/1963 
Rl as the definite members of subgroup I, while C/1882 
Rl, C/1945 XI, and C/1965 SI belonged to subgroup II. 
In addition, Marsden tabulated four and two additional 
comets, respectively, as possible members of the two cat- 
egories. A dent in this classification scheme was made by 
C/1970 Kl, whose orbit did not fit either subgroup. Be- 
cause its orbit was closer to subgroup II, Marsden (1989) 
classified it subsequently as subgroup Ila. Now comes 
C/2011 W3, whose perihelion distance fits subgroup I, 
but the angular elements are incompatible with any of 
the subgroups I, II, or Ila. This is worrisome, because 
Marsden (1967) regarded the longitude of the nodal line, 
one of the orbital angles, as the prime classifier for the 
subgroups. If one considers, by extension, C/2011 W3 to 
be a representative of a new subgroup III, one ends up 
with a total of eight bright sungrazers distributed into 
four categories, or two members per subgroup on the av- 
erage, not to mention that the orbit of C/1945 XI, one 
of the eight, is not really all that well fixed (sec Marsden 
1989 and our comment about Cunningham's computa- 
tions in Sec. 7). 

An overwhelming majority of SOHO, SMM, and SOL- 
WIND sungrazers belongs to subgroup I. This was shown 
to be the case by Marsden (1989) for the SMM and SOL- 
WIND comets and is illustrated for 1565 SOHO sun- 
grazers in a plot of the longitude of the ascending node 
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Figure 17. Comets C/2011 W3, C/1970 Kl, C/1965 SI, and C/1843 Dl in a plot of tlic longitude of the ascending node against the orbit 
inclination for 1565 SOHO Kreutz sungrazcrs from the period January 1996 to June 2010. C/2011 W3 is identified by the square, C/1970 
Kl by the spadesuit mark, C/1965 SI by the filled circle, while the location of C/1843 Dl coincides with the dense concentration of the 
SOHO sungrazers and is in the center of the open oversized circle. Comet C/2011 W3 appears to be orbitally related to a cluster of up to 
nearly 30 SOHO sungrazers with the ascending nodes between 310° and 335° and with the inclinations between 130° and 137°. 



against the orbit inclination in Fig. 17. However, nearly 
100 SOHO sungrazers, or about 6 percent of the total, 
make up a second branch, which in the range of ascend- 
ing node longitudes from ~320° to at least ~345° lies 
about 6° below an extended arm of the subgroup I set. 
And some 4° to 5° below this second branch there is yet 
another cluster of nearly 30 SOHO sungrazers with their 
nodal longitudes between 310° and 335°. 

For the sake of comparison, we have also included four 
bright sungrazers in Fig. 17. Comet C/1843 Dl, of sub- 
group I, is right in the dense core of the SOHO sungraz- 
ers, while C/1965 SI (of subgroup II), C/1970 Kl (of 
subgroup Ha), and C/2011 W3 are all located off the 
main branch. A surprising finding is that the second, 
much less populated but still well defined branch of the 
SOHO sungrazers belongs to subgroup Ha (not II, as has 
generally been assumed) . There is practically no concen- 
tration at all of the SOHO objects near the location of 
C/1965 SI. On the other hand, the agreement between 
the positions of C/2011 W3 and a third, very sparsely 
populated branch is rather obvious. Thus, when it comes 
to the SOHO sungrazers, the orientation of the orbit of 
C/2011 W3 in space is certainly not unique. And while 
the three branches essentially merge into a common re- 
lation in the plot of the longitude of the ascending node 
against the argument of perihelion in Fig. 18, C/2011 W3 
is again surrounded by a number of SOHO sungrazers. 

On the whole, the classification of the bright sungrazers 
into the subgroups is somewhat questionable. A sugges- 
tion for its eventual abandonment is supported by our 
finding (Sekanina & Chodas 2007) that even following a 
single fragmentation event at large heliocentric distance, 
two fragments of the same parent can end up in orbits 
formally belonging to different subgroups. In fact, the 



progenitor's splitting far from the Sun into two super- 
fragments with widely different orbits is a prerequisite 
behind the idea of introducing the concept of subgroups 
in the first place (Sekanina & Chodas 2004). 

The point of contention with C/2011 W3 is thus not 
whether it could derive from a precursor belonging to an- 
other subgroup, but, rather, that no plausible candidate 
for a sungrazing comet was ever recorded in the first half 
of the 14th century. Hasegawa & Nakano (2001) consid- 
ered only one possible sungrazer, X/1381 VI, during the 
entire 14th century. All three comets between 1282 and 
1368 in England's (2002) list have a relatively low sun- 
grazer ranking, even though the author does not entirely 
exclude the possibility for each of them, including one 
seen in Europe in October 1314, to be a member of the 
Kreutz system. 

On the other hand, only if the cumulative effect in the 
orbital period from a number of fragmentation episodes 
at large heliocentric distances could, in an extreme and 
quite unlikely run of events, amount to as much as ~50 
years, could X/1381 VI have been a parent to C/2011 
W3 in the sense that one of the products of its tidal- 
assisted breakup in 1381 would have become the first 
precursor to C/2011 W3. 

Of course, a large number of comets recorded in the 
late 13th century and during the 14th century were cat- 
aloged by Ho (1962), by Hasegawa (1980), and by Kronk 
(1999). However, most of them fail to fit a Kreutz sun- 
grazer because of the time of the year, the location in 
the sky, the apparent motion, or the appearance, while 
for the rest of them the reported information is too 
vague for the identification. The future modeling of the 
Kreutz-system evolution should proceed with these cir- 
cumstances taken into account. 
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Figure 18. Comets C/2011 W3, C/1970 Kl, C/1965 SI, and C/1843 Dl in a plot of the longitude of the ascending node against the 
argument of perihelion for 1565 SOHO Kreutz sungrazers from the period Jaimary 1996 to June 2010. The four bright comets are identified 
the same way as in Fig. 17. Although the three branches now essentially merge into one relation, the position of C/2011 W3 in the plot 
appears to be again closely matched by a modest number of SOHO sungrazers. 



14. CONCLUSIONS 

The issues and results presented in this investigation 
are summarized into the following conclusions: 

(1) Having survived the perihelion passage, comet 
C/2011 W3 was observed to undergo major morpholog- 
ical changes, culminating in the permanent loss of the 
nuclear condensation between 2011 December 19.4 and 
20.3 UT, some 4 days after perihelion. 

(2) This transformation was accompanied by a parallel 
development of a narrow, rectilinear dust streamer, ex- 
tending essentially along the tail's axis in the images of 
December 19 and rapidly evolving, in less than 24 hours, 
into a ribbon-like spine tail. 

(3) From the temporal variations in the projected ori- 
entation of its axis between December 19 and January 
18, we find that this spine tail was a synchronic fea- 
ture, whose brightest part consisted of submillimeter- 
sized dust grains released at velocities not exceeding 30 
m s~^ and whose origin was an outburst or a terminal 



(cataclysmic) fragmentation event (possibly even a rapid 
sequence of such very brief episodes), which peaked on 
December 17.6 ± 0.2 UT. 

(4) Accordingly, the comet's postperihelion activity, 
which lasted less than 48 hours and resulted in releas- 
ing up to an estimated 10^^ grams of dust (equivalent 
to a sphere 150-200 meters across), was the sole source 
of the spectacular dust tail, which was observed for 3 
months, until mid-March. 

(5) The delayed response of C/2011 W3 to the ex- 
tremely adverse environment in the close proximity of 
the Sun, which the nucleus' disintegration must have sig- 
nified, is of utmost importance. This overdue reaction 
speaks volumes about the cohesion of the nucleus and 
provides a clue to the probable nature of the breakup 
mechanism — as a product of a gradual heat pulse pen- 
etration to great depths of the nucleus and an entailing 
steady buildup of thermal stress throughout the interior. 

(6) Modeling of the heat-transfer process and the ther- 
mal stress distribution inside the nucleus, whose dimen- 
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sions diminish with time because of continuing erosion 
at the surface, shows that the rate of propagation of the 
heat pulse around perihelion was extremely rapid and 
highly asymmetric in that at greater depths beneath the 
surface both temperature and thermal stress were much 
higher at the same heliocentric distance after than be- 
fore perihelion and continued to grow far from the Sun, 
on the way to aphelion. 

(7) Based on the elastic constants for meteoritic, lunar, 
and some terrestrial analogs, the magnitude of thermal 
stress throughout much of the nucleus is estimated at 
more than 10 kPa shortly after perihelion, overcoming 
the object's expected cohesive strength (on the order of 
several kPa according to Lasue et al.'s aggregate model), 
riddling the entire interior with a dense network of clefts 
and cracks, thus causing eventually the nucleus' collapse 
into a cloud of debris. 

(8) The same mechanism may be responsible for the 
ubiquitous process of cascading fragmentation that sun- 
grazers are subjected to at all heliocentric distances, in- 
cluding very far from the Sun. 

(9) On the other hand, survival in the solar corona and 
the delayed response to the adverse conditions in its envi- 
ronment is inconsistent with the model of a strengthless 
or poorly cemented rubble-pile nucleiis made up essen- 
tially of unconsolidated debris. 

(10) If there were isolated reservoirs of volatile ices 
deep in the interior of the nucleus, major cracks caused 
by excess thermal stress would have opened them to the 
propagating heat wave, with the ensuing explosion con- 
tributing to the destruction of the nucleus; our model 
predicts that 1.6 days after perihelion, pockets of amor- 
phous water ice more than 40 meters beneath the surface 
would have been exposed to a temperature of ~130 K and 
thereby an exothermic reaction would be triggered, asso- 
ciated with amorphous ice's crystallization; only data on 
temporal variations in the comet's production of water 
vapor and/or other gas species could settle this issue. 

(11) The disappearance of the nuclear condensation 
deprived observers of obtaining postperihelion astrome- 
try needed for an accurate orbit computation by tradi- 
tional techniques, so that the orbital period of C/2011 
W3 remained essentially indeterminate even at the time 
of termination of observations. 

(12) We argue that the missing nucleus must have been 
located on the synchrone defined by the axis of the spine 
tail whose orientation and sunward tip have been mea- 
sured; the missing nucleus also must have been situated 
on the line of forced orbital-period variation, computed 
from the orbital solutions based on high-quality preper- 
ihelion astrometry from the ground. 

(13) We succeed in deriving the astrometric positions 
of this missing nucleus as the coordinates of the points of 
intersection of the spine tail's synchrone and the line of 
orbital-period variation; a high-quality orbital solution is 
obtained by linking the preperihelion observations with 
these derived positions of the missing nucleus. 

(14) The resulting orbit gives 698 ± 2 years for the os- 
culating orbital period and confirms that C/2011 W3 is 
the first major member of the expected new, 21st cen- 
tury cluster of bright Kreutz-system sungrazers, whose 
existence was predicted by these authors in 2007. 

(15) From the spine-tail evolution we determine that 
the measured sunward tip receded antisunward from the 



computed position of the missing nucleus, the distance 
increasing from 1' on December 23 to more than 10' in 
mid- January and that this terminus was populated by 
dust particles 1-2 mm in diameter. 

(16) The bizarre appearance of the comet's dust tail in 
images taken with the coronagraphs on board the SOHO 
and STEREO spacecraft only hours after perihelion is 
modeled, using the standard dynamical method, in terms 
of the ratio j3 of the acceleration by solar radiation pres- 
sure to that by the solar gravitational attraction; the tail 
brightness was affected by forward scattering of sunlight 
by dust to a considerable degree. 

(17) Modeling of the SOHO's C3 preperihelion images 
of the dust tail shows a population of predominantly 
microscopic particles with a broad range of sizes and 
composition, described both by /3 > 1 (strongly absorb- 
ing, such as carbon-rich, grains moving in convex hyper- 
bolic orbits) and by /3 < 1 (essentially dielectric, perhaps 
mostly silicate, grains in concave hyperbolic orbits); the 
ones that were still imaged after perihelion had /? ~ 0.6, 
probably submicron-sized silicate particles. 

(18) By contrast, nearly all microscopic dust released 
before perihelion into convex hyperbolic orbits moved al- 
most immediately away from the Sun and was not seen 
in the postperihelion images of the comet. 

(19) The disconnection of the comet's head from the 
dust tail released before perihelion, strikingly depicted 
in the postperihelion images taken with the SOHO's 
C2 coronagraph and with the coronagraphs on board 
STEREO- A and STEREO-B, is interpreted as a result 
of vigorous sublimation of submicron-sized dust at helio- 
centric distances smaller than about 1.8 Rq, primarily 
particles with (3 < 0.6 released between 0.1 day before 
and 0.1 day after perihelion. 

(20) Prom the integrated effect of dust particle sub- 
limation, we establish the relationship between the 
sublimation-pressure constant A and the heat of subli- 
mation L and find that it is consistent with the L{A) 
relation for Mg-rich olivine and that therefore it is likely 
that the dust sublimating very close to the Sun is domi- 
nated by olivine-based silicates, a conclusion that is also 
supported dynamically by the magnitude of the acceler- 
ation parameter /3. 

(21) The place of C/2011 W3 in the hierarchy of 
the Kreutz system of sungrazing comets is a matter of 
some speculation, but one possible evolutionary path 
is charted, based on our orbit integration back to the 
comet's previous return to perihelion, which nominally 
occurred in 1329. 

(22) In the context of Scenario B of the Kreutz sys- 
tem's hierarchy (from our 2007 paper) we follow the tidal- 
assisted breakup of the progenitor (or superfragment) — 
the comet of 467 and the expected chain of nontidal 
fragmentation events into a number of parent fragments 
(including X/1106 CI); a subsequent tidal-assisted split- 
ting of one such parent fragment in the 14th century into 
another generation of fragments, including the first pre- 
cursor of C/2011 W3; and, finally, another bout of non- 
tidal cascading fragmentation with the birth of comet 
C/2011 W3 discrete object of its own. 

(23) In this scenario, C/2011 W3 is indirectly related 
to X/1106 CI and C/1843 Dl, and in the coming years 
and decades it may be followed by more equally bright 
or brighter sungrazing comets in similar orbits. 
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(24) Comet C/2011 W3 presents a complication for the 
classification of the bright sungrazers into subgroups: its 
perihelion distance is typical for subgroup I, but its an- 
gular elements do not fit any of the subgroups I, II, or 
Ila; thus, among the bright members of the Kreutz sys- 
tem, the orbit of C/2011 W3 is unique, although we are 
not ready to call it a representative of a new subgroup. 

(25) In a plot of the longitude of the ascending node 
against the inclination, C/2011 W3 is associated with 
the least populated of the three branches into which the 
SOHO sungrazers discriminate; from this standpoint, its 
orbit is not unique. 

(26) Because the orbits of the sungrazers (including 
the bright objects) can be transformed by the cascading 
fragmentation process from one subgroup into another 
subgroup during one revolution about the Sun, the clas- 
sification of the bright sungrazers into subgroups may no 
longer be appropriate. 

(27) Future modeling of the hierarchy and evolution of 
the Kreutz system of sungrazers should account for the 
existence of C/2011 W3 and its apparent relationship 
with the minor branch of SOHO sungrazers; at a 
minimum, orbital computations should verify the degree 
of plausibility of a fragmentation path of the type 
displayed in Fig. 16. 
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